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Abstract 


A novel method is developed for producing crystalline 
silicon pn homojunctions, using doped amorphous silicon 
as source for the dopant diffusion. A highly phosphorus 
doped thin layer of hydrogenated amorphous silicon was 
pre-deposited by plasma enhanced chemical vapour de- 
position on a p-type doped c-Si wafer. Samples were then 
heated to high temperatures for the P to diffuse from the 
amorphous to the crystalline silicon, creating the pn ho- 
mojunction. A J(V) characteristic curve under AM1 equi- 
valent conditions was obtained, reading an efficiency of 
5.2 %. 


Work Overview 


Crystalline silicon (c-Si) wafer pn homojunctions are con- 
ventionally produced by a two-step process - i) the de- 
position of a dopant dose and ii) the dopant drive-in. The 
first step is normally made by one of two methods, the 
dopant pre-deposition at high temperature or the ionic 
implantation. In the dopant pre-deposition method, a thin 
layer of a highly concentration dopant compound origina- 
ted from a solid, liquid or gas source, is deposited on the 
c-Si by heating the wafer to high temperatures (>800 ºC), 
to increase the dopant solid solubility in silicon. The re- 
latively low production cost achieved using this method 
makes it the most used by the c-Si photovoltaic industry 
at the moment. On the other hand, the ionic implantation 
involves injecting the dopant directly into the wafer, whi- 
ch requires relatively expensive equipment, but also allo- 
ws for enhanced quality junctions and good process con- 
trol. In this work, an alternative method has been found 
for producing crystalline silicon homojunctions. 


Phosphorus doped hydrogenated amorphous silicon (a- 
-Si:H) was pre-deposited by Plasma Enhanced Chemical 
Vapour Deposition (PECVD) on a boron doped c-Si Cz 
<100> wafer (p = 1.6 Q.cm). The PECVD deposition was 
done using a constant total gas flow of 25 sccm compo- 
sed mainly of silane and phosphine (0.06 % of the total 
flow) during 1 min at RF power Ppr = 62 mW/cm? and 
substrate temperature T = 350º. The resulting film has a 
thickness of 58 nm, measured by profilometer, and hydro- 
gen content between 15 and 17 at.%, as determined in a 
previous study [1]. 
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The second step was to introduce the sample in a diffu- 
sion furnace and heat it at 5 “C/min up to 350 “C, from 
there up to 550 “C at 0.5 “C/min and then keep the sam- 
ple at this temperature for 1 h. This process allows for a 
smooth dehydrogenation of the a-Si:H sample, avoiding 
damage in the film, as studied in the past [1,2,3,4,5,6]. 
Finally, the sample was heated from 550 “Cup to 900ºCat 
5 “C/min and kept at the latter temperature for 1 h, while 
injecting a constant inert gas flow of N; at 1 1/s, in order 
to avoid undesired oxidation of the sample. At 900ºC, the 
phosphorus contained in the amorphous silicon thin film 
diffuses towards the c-Si wafer, consequently creating a 
pn homojunction. Dehydrogenating smoothly in the pre- 
vious step is important because a fast dehydrogenation 
would cause craters in the film, and posteriorly there 
would be no dopant diffusion in those regions, so that 
depositing a metallic contact above the junction would 
short circuit it. The full process, in three steps, is schema- 
tically summed up in figure 1. 


1. Pre-deposition of a-Si:H(n) 
by PECVD 


2. Dehydrogenation at 
350-550 “C 


3. Phosphorus diffusion from 
the a-Si to the c-Si at 900 C 


Figure 1.Schematic description of the doping of c-Si by 
pre-deposition of a-Si:H(n) 


In order to measure the photovoltaic J(V) characteristic curve of the obtained junction, a 3 mm disk of around 130 nm 
thickness of aluminium was deposited by resistive evaporation on top of the sample and then annealed during 30 min at 
400 ºC, in order to better collect current. Additionally, 140 nm of In>03 which is a transparent conductive oxide (TCO) was 
deposited by plasma enhanced reactive thermal evaporation, to further improve the current collection by the pn junction. 
The TCO surface is a 0.36 cm? square and covers completely the aluminium contact, so that the total area covered by the 
TCO and not by the aluminium is 0.29 cm?. The J(V) curve in the dark and under AM 1 equivalent conditions is plotted in 
figure 2, reading an efficiency of 5.2 %. 


J (mA/cm?) —— in the dark 
60 


—— under AM 1 equivalent 
conditions 


Efficiency : 5.2 % 

Fill Factor : 28 % 

Voc :0.54V 
Jsc : 34 mA/em 
R parallel : 20 O.cm? 
Rseries :130.cm? 


Figure 2. J(V) characteristic curve of the obtained pn junction in the dark and under AM 1 equivalent conditions 


The obtained curve shows that this pn homojunction production method, alternative to the conventional pre-deposition 
and ionic implantation, allows for the production of solar cells. Furthermore, higher efficiencies should be achieved by 
improving the series resistance, which here seems to be masking an ideal curve. This may be due to the a-Si layer, which 
has only partially crystallized during the drive-in phase, or alternatively to the contact resistance between the n-type 
doped layer and the TCO. The parallel resistance, which was calculated at | = O mA/cm2, seems to be masked by this high 
series resistance. As a result of these, an optimization of the method and additional work on getting better contacts 
should improve the efficiency of such a solar cell, eventually making it an interesting option for industrial purposes. 
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Abstract 


Thermal electricity generation has great importance on the security 
of electricity grids and also in energy costs control, 

n Portugal, last years, coal plants supplied an average of 22,3 % of 
consumption. 
But burning coal raises environmental problems, such as air emissions, 
which requires the introduction of specific technologies to minimize 
its effects. 
This includes dust removal, desulphurisation and denitrification. 
Products resulting of these equipments operation, fly ash, bottom ash 
and FGD gypsum, should be disposed in large landfills, if nothing was 
done, 
The aim of this document is to show that, with the right management, 
coal combustion products can no longer be considered as industrial 
waste but, instead, they are raw materials, with commercial value. 
Keywords: Electricity generation, thermal power plants, coal, 
environment, coal combustion products, fly ash, bottom ash, gypsum. 


Resumo 


A produção de electricidade com origem térmica tem uma importância 
muito grande na segurança das redes eléctricas e também no controlo 
dos custos da energia. 
Em Portugal, nos últimos anos, as centrais a carvão satisfizeram, em 
média, 22,3 % do consumo, 
Mas a queima do carvão levanta problemas ambientais, a nível 


de emissões atmosféricas, cuja minimização exige a introdução 
de tecnologias como o despoeiramento, a dessulfuração ou a 


desnitrificação. Do funcionamento destes equipamentos resultam 
resíduos industriais como as cinzas, volantes e de fundo, e o gesso, 
A sua existência levaria, se nada fosse feito, à sua deposição em 
grandes aterros. 


o presente trabalho demonstra-se que, através duma gestão 


dequada, os produtos da combustão do carvão podem ter utilizações 
ue os retiram da sua classificação de resíduos. 


a 
q 
Palavras-chave: Produção de electricidade, centrais térmicas, carvão, 
ambiente, cinzas, gesso. 
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1+ INTRODUÇÃO 


Na produção de electricidade as tecnologias associadas ao carvão 
desempenham um papel relevante na satisfação dos consumos. 

O carvão é a ºº* fonte de abastecimento de energia primária a nível 
mundial, com o peso de 30,3% em 2011, só superado pelo petróleo 
com 33,1%. 
Daí que, sendo reconhecido o impacto ambiental da utilização das 
energias fósseis, no caso do carvão é de extrema importância fazer a 
gestão dos produtos resultantes da sua queima, 


E se nos gases essa gestão atravessa uma fase de inovação e 
investigação, no caso dos produtos sólidos a solução elementar é a 
sua deposição em aterro. 


Este trabalho pretende demonstrar a importância da actividade de 
utilização dos produtos da combustão do carvão (PCC), apresentando 
as soluções adoptadas. 


o seu desenvolvimento segue-se a cadeia de produção; 
— O carvão - a electricidade - os resíduos - a utilização 


Está presente a regra dos 3R, que geralmente se invoca para lidar com 


o problema dos resíduos: 

Reduzir, com as soluções correspondentes a um melhor desempenho; 
Reutilizar, no caso da recolha das cinzas de fundo e das tecnologias 
de beneficiação; 

Reciclar, aproveitando os produtos, que se tornam matérias-primas 
de outros processos de produção. 


2 * O CARVÃO 


As reservas provadas de carvão em 2011 [1] eram suficientes para a 
satisfação do consumo durante 112 anos, sendo assim o combustível 
fóssil mais abundante, pelo rácio reserva/consumo. 

Este rácio é de 54,2 anos para o petróleo e 63,6 anos para o gás 
natural, 

Para o conjunto da União Europeia tem-se vindo a verificar ao longo 
do tempo uma tendência de descida do consumo dos diversos tipos 
de carvão, como se podem ver pela consulta das bases de dados da 
Eurostat. 


Esta redução é contrariada nos 2 últimos anos, de 2010e 2011, onde 
se verifica um acréscimo de consumo. 

Tal facto é atribuído à competitividade de preços com os outros 
combustíveis e ao incremento do uso do carvão na produção de 
electricidade em países como a Alemanha e a Espanha. 


2.1 * O CONSUMO DE CARVÃO EM PORTUGAL 


Ao longo dos últimos anos é visível o decréscimo do petróleo, cada 
vez menos usado na produção de electricidade e a ascensão das 
renováveis, impulsionadas pela energia eólica [Fig. 1]. O gás natural 
cresce principalmente no sector doméstico, remetendo-se assim o 
carvão para uma situação de apoio, o que lhe valeu a recuperação em 
2011 eem 2012, resultante de serem anos secos, 


Quantidades (ktep) 
o 
o 
[=] 
[=] 


2000 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 


Anos 


— petróleo — Carvão — Gás natural — Renováveis 


Figura 1 - Consumo de energia primária 
(Fonte: DGEG, Balanços Energéticos) 


No que diz respeito à utilização, existem 3 grandes sectores onde o 
carvão é o combustível preferencial: 


- O siderúrgico e metalúrgico, que em Portugal foi perdendo peso 
como consumidor desta forma de energia, designadamente com a 
reconversão da Siderurgia Nacional em que o alto-forno foi substituído 
por um forno de arco voltaico. 


- À indústria cimenteira, mas que face aos preços vantajosos do 
petcoke, passou a dar preferência a este combustível. 


- A produção de electricidade, que representa assim a quase totalidade 
dos consumos nacionais. 

As centrais térmicas a carvão 

A partir de informação pública disponibilizada pelas empresas, 


descrevem-se as principais características das centrais térmicas a 
carvão, actualmente em serviço em Portugal, 


Figura 2 - Central de Sines 
(Foto: EDP) 


Características gerais [2] 

Empresa: EDP Gestão da Produção de Energia, SA 
Concelho: Sines 

Entrada em serviço: 1985, 1º grupo; 1989, último grupo 
Potência à emissão: 1180 MW 

Número de grupos: 4 


Figura 3 - Central do Pego 
(Foto: Tejo Energia) 


Características gerais 

Empresa: Tejo Energia, SA 

Concelho: Abrantes 

Entrada em serviço industrial: 1993, 1º grupo; 1995, 2º grupo 
Potência à emissão: 5/6 Mw 

Número de grupos: 2 
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2.2 « O CARVÃO EM PORTUGAL 


Na consulta das estatísticas da DGEG [Tabela 1] é possível ver o 
carvão como recurso energético endógeno, até ao ano de 1994, 
altura em que cessou a produção das minas de S, Pedro da Cova e do 
Pejão, cujo carvão abastecia a Central Térmica da Tapada do Outeiro. 
Até ao seu encerramento, em 2004, esta Central passou a consumir 
exclusivamente fuelóleo, 


Tabela 1 - Produção de carvão nacional (ktep) 


Ano: 1990 1991 1992 1993 1994 1995 
Carvão 115 111 91 81 60 0 


O facto de existirem reservas deste combustível fez com que ele 
desempenhasse uma função importante no abastecimento energético 
do país, com um papel relevante na história do sector. 


3 « À ELECTRICIDADE 


Num sistema concorrencial, como o Português, actualmente, a alocação 
da energia eléctrica, mesmo que os grupos não obedeçam a um 
despacho económico, só despacho técnico para garantir a segurança 


ra 


da rede, não é muito diferente do despacho a custo marginal, 


Com efeito, entrando em serviço os grupos que “casaram” a oferta e a 
procura, significa que entram primeiro os grupos mais baratos, sendo a 


ponta assegurada pelos mais caros, tal como no despacho económico 


a custo marginal, 


3.1 O PAPEL DAS CENTRAIS A CARVÃO 


Com a penetração dos produtores em regime especial, sem actuarem 
no mercado e com a sobrecapacidade do sistema, as centrais a carvão, 
desempenham tendencialmente o papel de fecho do diagrama de 
cargas: predominando nas pontas, entram depois das PRE, fios 
de água, ciclo combinado (antes do “take or pay” contratual do gás 
natural) e importação. 


Claro que esta utilização tende a aumentar em anos secos, como 2008, 


2009 ou mesmo em 2012, mas com a hidraulicidade a aumentar o seu 


funcionamento diminui [Tabela 2]. 


Tabela 2 - Satisfação dos consumos pelo carvão (Fonte: REN) 


Carvão 


Observa-se que o ano de 2008 foi afectado pela indisponibilidade 
associada à instalação dos equipamentos de dessulfuração. 

No futuro, as centrais a carvão poderão ter um funcionamento mais 
estável, aumentando o peso das hidroeléctricas com bombagem ou, 
se a isso ndo corresponder uma desaceleração de energias como a 


eólica, podem definitivamente ficar remetidas para a inactividade. 


Claro que aqui também entra o valor relativo dos custos de produção a 


gás natural e a carvão, onde tendencialmente este último combustível 
tem levado vantagem. 


3.1 « O CARVÃO LIMPO 


Em Portugal, quando por opção de política energética para a 
produção de electricidade, foram construídas as centrais térmicas a 
carvão de Sines e Pego, a minimização do impacte ambiental foi uma 
preocupação patente nos projectos e nas actualizações que sofreram 
posteriormente com a introdução das tecnologias de dessulfuração e 
de desnitrificação [Fig. 4]. 


Amónia Chaminé 
NH, 


Desnitrificação 


É 

E 

& 

Calcário) m 

—— E 

TT I 
Cinzas de fundo Cinzas volantes Gesso 


Figura 4 - Gerador de vapor de uma central a carvão 
(Fonte: ECOBA, adaptado do inglês) 


A tecnologia utilizada nestes Grupos Geradores de Vapor leva à 
produção de três produtos resultantes da combustão do carvão e do 
processo de dessulfuração: 
- As cinzas volantes, recolhidas nos precipitadores electrostáticos; 

- As cinzas de fundo obtidas por gravidade nos cinzeiros das caldeiras; 


- O gesso resultante da captura do óxido de enxofre nos gases. 
Foi prevista uma área de deposição junto às centrais para todos 


estes produtos. No entanto sempre existiu a preocupação de que a 


colocação em aterro fosse o último recurso, por motivos económicos 
e ambientais. 


são pa batcmet de nsumo | 2005 2506 [2507 [ 200 [ 2008 [ amio [am 
Ses) | 550/ S64/ G00[ 65%[ su0/ som| com 


Pego Pego(cwh) 


[ATo[ 4576] S015[ a56/ 3075 voss/ 220 


a — [Ns avo[is[ os rae assar 
[mecnitom | 9750/ 19774] SO059[ SO505[ A507A/ sej99/ 50503 
Perdonbo] | m8[ =| =| ms| =5| v6| ot 
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Quando a redução das emissões dos gases com efeito de estufa, 
alegadamente responsáveis pelo aquecimento global do planeta, se 


nFA “BA Bs «FEC =SDA «FGD 


tornou um dos objectivos fundamentais das políticas ambientais, o 
C0> ficou no centro das atenções [3]. A sua captura e armazenamento 
é uma área onde têm sido concentrados esforços de investigação. As 


principais tecnologias testadas são: post-combustão, pré-combustão 


e oxi-combustão. 


CCP production in milion tormes 


4 e AS UTILIZAÇÕES 
Os produtos resultantes da combustão de carvão em centrais PPP 
sctri 3 icã » ] PLS) 1) a! o 
termoeléctricas exigem uma gestão que evite a sua deposição. PPP SDS SDS 
A cooperação entre produtores e indústria, cimenteira e de fabrico FA = iy as; BA - Donom ash; ES = bober siag: FC - FEC asi. SOA - SDA proguct; FGD - FGD gypsam 
de betão, tem conduzido a soluções sustentadas e vantajosas para 
ambas as partes. Figura 5 - Produção dos PCC 
Do lado da produção, é impossível desligar estes produtos das (Fonte: ECOBA) 


condicionantes a que ao sector energético está sujeito. 


Inversamente, sendo a sua utilização feita quase exclusivamente no 
sector da construção civil, também do lado da procura, a conjuntura 


económica é determinante. Utilisation in 
Construction 
Industry and 


Underground 


4.1 « SITUAÇÃO NA EUROPA Mining 


Temporary 
Stockpile, 2.2% 


N Disposal, 6.4% 
Em 2010, cerca de 49 milhões de toneladas de PCC foram produzidos 
na Europa (UE15). 

A produção em todos os estados membros da Comunidade Europeia 


(CE) é estimada em cerca de 100 milhões de toneladas. 
Restoration of 

Open Cast Mines, 
Quarries and Pits 


Os valores globais da Europa dos 15, em 2010, encontram-se 
compilados em [4] 

Interessa analisar alguns aspectos e tendências mais relevantes [Fig. 
5e 6]: Figura 6 -Utilização dos PCC 

(Fonte: ECOBA) 

Desde 2003 que a produção de PCC tem-se vindo a reduzir fruto 

da menor utilização das centrais a carvão por motivos de menor 
actividade económica (os consumos de electricidade têm-se reduzido), 

mas também ambientais e da introdução de energias renováveis 
No capítulo da utilização verifica-se a existência de escoamentos 
estáveis sendo que o recurso à deposição representou só 6,4% da 
produção. Os produtos que são depositados em aterro dificilmente 
encontrariam mercado por falta de enquadramento normativo. 


As grandes utilizações são correspondentes às cinzas, cujos principais 
destinos são a construção e a recuperação mineira. 


No que diz respeito ao gesso [Fig. 7] é evidente o resultado das 


FGD gypsum production 
in million tonnes 


tecnologias de dessulfuração, nos diferentes países da Europa, através 


de actualizações nas centrais já existentes ou no equipamento das 


1993 1905 1997 1999 2001 2003 2005 2007 2009 


novas instalações. 


A produção sobe de um modo muito significativo a partir do ano 2000, Figura 7 - Produção de gesso, EU15, 1993-2009 
quer para as centrais de carvão betuminoso, quer para as de queima (Fonte: ECOBA) 
de lenhites. 
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4.2 e SITUAÇÃO EM PORTUGAL 


A partir dos dados publicados e por estimativa dos restantes, 


consegue-se caracterizar a situação nacional dos PCC, até 2011 
[Tabela 3]. 


CE ESA E E 
Quanto à utilização, relativamente às cinzas 


7536] 147345 | 241156] 135221| 170 132 volantes, o escoamento é feito, 
TOTAL (t) | 456 202 | 577941| 756861| 410622 | 516 102) em média, por [Tabela 4]; 
Média 2005-2011 Ya 
Tabela 3 - Produção de PCC em Portugal 


As cinzas de fundo só esporadicamente têm escoamento, como Adição de betão 
adição de betão ou constituinte minoritário do cimento. Enquanto 
a construção rodoviária não tiver valor expressivo, é feita a sua 


O gesso tem vindo a adquirir mercado, principalmente no fabrico de 
placas. O facto de a produção ter vindo a encontrar esta utilização, Tabela 4 - Utilização de cinzas volantes 


quer no mercado nacional, quer para exportação, leva a que as 


quantidades depositadas o sejam como armazenagem temporária. Geotecnia [6-9] 
Tendencialmente, os respectivos aterros vão ficando sem material 
depositado. Na construção de estradas, as cinzas volantes e de fundo são 
particularmente adequadas como material de enchimento nas bases 
4.3 « PRINCIPAIS APLICAÇÕES e sub-bases dos pavimentos 
. Numa fase em que todos os factores de redução de custos se 
Indústria cimenteira e de betões tornam importantes, o recurso a materiais alternativos, de qualidade 


comprovada, é um factor que aumenta a competitividade para aquelas 


Este enquadramento deixa antever que as cinzas volantes se empresas que se dispõem a assumir um papel na inovação. 


encontram ligadas a um conjunto de importantes obras públicas que 


constituem referências do sector em Portugal, Alguns exemplos [5]: 
Gesso 


- Blocos do molhe do porto de Sines, edifício sede da Caixa Geral de 
Depósitos, via do Infante, metro do Porto; Com a instalação em Portugal de sistemas de dessulfuração, em 


Edo 2007/2008, os operadores das centrais térmicas necessitaram de 
- Estádios do Campeonato Europeu de Futebol, 
encontrar um escoamento para o gesso. Em 2011 foram produzidas 


- Ponte Vasco da Gama, ponte da Lezíria, barragens dos novos cerca de 170 000 +. 


aproveitamentos hidroeléciricos Como factores favoráveis surgiram também 2 realidades, que 


O interesse comum a estas utilizações reside nas propriedades prevaleceram: 


especiais que o betão apresenta e que são nomeadamente: «Oigêsso produzido revelbu-se de excelente qualidade 


- Progressão do endurecimento durante um período longo; - O mercado nacional das placas cartonadas e das placas de gesso 


- Bom comportamento em ambientes agressivos; aminadas, estava em desenvolvimento. 


- Elevada compacidade e impermeabilidade; ar e = 
as na prática a solução passa pelas aplicações mais comuns como 


- Ganhos na trabalhidade e nas condições de aplicação, seja na grandes consumidoras: 


bombagem, seja na compactação; 
- Regulador de presa no cimento 


- Menor necessidade de água na amassadura, A a 
- Gesso projectado para acabamento de interiores 


- Placas de cartão-gesso 


- Blocos de gesso 
- Argamassa autoniveladora de pavimentos 
- Aterros e enchimentos 


5 « CONCLUSÕES 


Com base nas directivas e nos regulamentos que afectam a produção 
de energia nos estados-membros da União Europeia, - e, portanto, 
também a produção de PCC - a indústria de energia terá de desenvolver 
todos os esforços para fornecer sempre produtos de boa qualidade 
para o mercado da construção. 

A gestão desses produtos de acordo com os padrões de certificação, é 
uma garantia da qualidade exigida, 


Citam-se de seguida 3 vertentes desta gestão que ilustram factores 
de contributo decisivo para o interesse económico e ambiental, seja 
pela perspectiva da economia de recursos naturais e energéticos, seja 
pela redução de emissões atmosféricas: 


A normalização 


Integrando processos de fabrico, como o cimento e o betão, que estão 
sujeitos a regulamentos, normas, especificações muito detalhadas, as 
cinzas têm por sua vez de ser também certificadas. 

No caso do gesso, o que é importante é o cumprimento dos requisitos 
da Eurogypsum. 
Assim os geradores de vapor e os sistemas de controlo de processo 
foram concebidos para permitirem que as cinzas volantes tivessem de 
acordo com EN 450 de “Cinzas volantes para betão”, 

Desde 2006 que em Portugal as cinzas volantes, através de 


certificação, possuem a marca CE, o que permite que em caso de 


utilização transfronteiriça, no espaço europeu, a certificação seja 
válida. Isto é, as cinzas Portuguesas podem ser vendidas em qualquer 


país Europeu. 


Mesmo as cinzas de fundo extraídas em cinzeiro seco, têm todas as 


condições para ser utilizadas para os mesmos fins das cinzas volantes, 


se misturadas com elas, até ao máximo de 10% em massa. 


A classificação dos PCC 


Actualmente existem mecanismos legais que permitem que, apesar 
da classificação de resíduo, estes materiais possam ser equiparados 
a subprodutos. 

Foi o que aconteceu com as cinzas volantes no betão e no cimento, 
há alguns anos, 

É a situação do gesso para placas cartonadas e placas laminadas. 


São claros os critérios a que devem obedecer e estão mencionados no 
artº 5 da Directiva Quadro dos Resíduos [10]: 

- À garantia da sua utilização futura; 

- À possibilidade de utilização directa, sem processamento posterior; 
- Ser parte integrante de um processo de produção; 

- O seu uso estar de acordo com a legislação ambiental e de saúde, 
sem impactes adversos, 

Depois do registo no REACH, fazia sentido a desclassificação destas 


substâncias como resíduos, passando a subprodutos invocando como 


argumento os trabalhos e os testes necessários para aquele processo. 


A cooperação 


O relacionamento entre as entidades que possuem o mesmo tipo de 
problemas e assim podem partilhar soluções é o fundamento de uma 
política sempre adoptada pelos produtores europeus de electricidade. 
A participação activa na Associação Europeia dos Produtos da 
Combustão do Carvão, ECOBA, de que a EDP é membro fundador, 


é uma manifestação desta atitude e tem levado à obtenção de 
vantagens mútuas. 


Por outro lado, a experiência adquirida ao longo dos anos na gestão 


te tipo de materiais permite uma análise dos factores que 


am a sua utilização e a identificação das medidas que se podem 


u 
implementar para reverter a situação, 
E uma análise que passa por questões operacionais, legais e 


económicas, pois o aspecto ambiental apresenta vantagens 


aparentemente inquestionáveis. 
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Resumo 


O estudo de buracos negros é uma área de investigação da 
física teórica de intensa atividade. Neste artigo analisar-se-á 
um método para o estudo da entropia destes corpos baseado 
no conceito de camadas finas de matéria. 


1 + INTRODUÇÃO: BURACOS NEGROS 


Em 1916, Einstein generalizou a gravidade de Newton para um con- 
texto relativístico, em que a velocidade da luz não só existe como 
limite inatingível de velocidade para qualquer corpo como também 
é localmente constante para todos os observadores. Nesta teoria, 
o tempo e espaço deixam de ser absolutos, entrelaçando-se antes 
numa entidade designada de espaço-tempo. A atracção gravitacio- 
nal é então explicada como sendo uma consequência da curvatura do 
espaço-tempo gerada por corpos massivos. A força gravítica, que an- 


+ 


eriormente era vista como uma acção à distância, passou a ser com- 


+ 


preendida como um simples movimento natural dos corpos no espaço- 


empo curvo, em trajectórias designadas de geodésicas. 


O desenvolvimento desta nova teoria veio acompanhado com um novo 


ormalismo matemático. O papel do potencial gravítico passou a ser 
desempenhado pela métrica, uma noção com origens na geometria 
diferencial, que contém essencialmente todas as propriedades que 
caracterizam a geometria do espaço-tempo. Os objectos matemáticos 


de interesse passaram a ser os tensores, uma generalização do con- 


ceito de vector, e a lei da gravitação de Newton foi substituída pelas 
equações de Einstein 


Gab =BnTap (1) 


que relacionam o tensor de Einstein Gap - que é função da curvatura 
do espaço-tempo e envolve derivadas segundas da métrica em re- 
lação ao espaço-tempo - com o tensor Tap de energia-momento da 
matéria presente. Note-se que na equação (1) foram usadas unidades 
naturais em que a velocidade da luz ce a constante de gravitação uni- 
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versal G foram postas iguais a 1, Se assim não fosse, a constante de 
proporcionalidade seria 8xG/c? em vez de 8m. 


As equações de Einstein são complexas o suficientemente para pos- 
suírem soluções exactas de elevado interesse físico. Uma delas é a 
solução de Schwarzschild, encontrada em 1916, dada pelo elemento 


de linha 


-1 
ds2=-(1-Dar+ (1) dr? +rê(do? +sin?0 dg?) (2) 


r 


onde m é a massa geométrica, dada em termos da massa usual M por 
m=MG/c2. A coordenada t é a coordenada temporal e (1,8,4) são as 


E 


oordenadas esféricas usuais. A expressão (2) descreve a curvatura de 
um espaço-tempo esfericamente simétrico no vácuo, e apresenta um 
comportamento peculiar para valores de coordenadas radiais menores 


U iguais a R.=2m, chamado raio de Schwarzschild, raio gravitacional 


o) 
ou raio de horizonte de eventos. Nesse próprio raio a métrica parece 


singular mas não o é e, para raios menores, as componentes da mé- 


trica temporal e radial trocam de sinal, invertendo consequentemente 


o comportamento dessas duas coordenadas. Acontece que essa apa- 
rente singularidade no raio de horizonte deve-se puramente à escolha 


de coordenadas, sendo prontamente removida se for feita uma es- 


colha mais adequada. No entanto, é de salientar que esse raio cor- 


responde fisicamente a um ponto sem retorno, já que se uma partícula 
entrar na esfera de Schwarzschild - esfera de raio Rç - seguirá inexo- 
ravelmente para a singularidade central, essa uma singularidade real 


nde a curvatura é infinita. À região contida dentro do horizonte de 


o) 
eventos dá-se o nome de buraco negro, 
o) 


Buracos negros são o resultado do colapso gravitacional de estrelas 


u de grandes aglomerados de matéria, podendo ter massas da ordem 
de 1 MO, onde a notação MO corresponde a uma massa solar, até 


massas da ordem de 10º MG, aos quais se dá o nome de buracos 
negros supermassivos. 


Este último tipo de buracos negros encontram-se no centro de qua- 
sares e de galáxias. A nossa própria galáxia tem um buraco negro cen- 


ral com uma massa da ordem 10º MO. 


O estudo da interacção de buracos negros com matéria e radiação 
elecromagnética revelou a importância destes objectos dentro da físi- 
ca, nomeadamente quando efeitos quânticos são levados em con 


(ab) 


na vizinhança do horizonte de eventos de um buraco negro. Nes 


(ab) 


situação, uma multitude de novos aspectos notáveis aparecem, en- 
tre um deles a radiação de Hawking. Este tipo de radiação pode ser 
interpretada como emergindo do princípio da incerteza de Heisen- 
berg que afirma que num intervalo de tempo suficientemente curto 


a incerteza na energia pode ser grande o suficiente para aparecer um 


4 


ti-partícula, desde que o par se aniquile nesse 


par de partícula e an 


intervalo de tempo. Se isto acontecer perto do horizonte de eventos, 


uma das partículas do par criado pode ser capturada pelo buraco negro 


enquanto que a outra escapa. Um observador exterior ao horizonte 


f 


de eventos percepciona este efeito como a emissão de uma partícula 
por parte do buraco negro com a consequente perda de energia por 
este, A radiação de Hawking atribui assim uma temperatura a buracos 


negros, chamada temperatura de Hawking, mostrando que estes se 


comportam como sistemas termodinâmicos, de facto corpos negros, 
aos quais está naturalmente associada uma entropia. Essa entropia é 


dada pela fórmula de Bekenstein-Hawking [3,4], nomeadamente 
Son = (3) 


onde A é a área do horizonte de eventos e onde se usam as uni- 
dades naturais c=G=h=1, em que h=h/27 e onde h é a constante de 
Planck. Se recuperarmos estas constantes universais, a fórmula fica 


14 
San = 4 ap Planck (4) 


onde A pl é a área de Planck dada em termos das unidades fundamen- 
tais por Gh/c3 e tem o valor de 2.61x10-76 cm? e Spunck” Kg Onde kg 
é a constante de Boltzmann. Esta área é formada pelo quadrado de 
comprimento de Planck, considerado o menor comprimento com sig- 
nificado físico. 


A expressão (3) é especialmente importante pelo facto de surgir de 
uma combinação de relatividade geral com mecânica quântica já que 
entram as contantes universais correspondentes a cada uma dessas 
teorias, constituindo assim um resultado que uma teoria de gravita- 
ção quântica deverá conter no limite semiclássico. Desta forma, é de 


esperar que um estudo mais profundo da entropia de buracos negros 


possa indicar novas pistas para o desenvolvimento de uma teoria de 
gravitação quântica. 


Para uma análise mais detalhada de buracos negros sugerem-se as 
referências [1, 2], para a temperatura e entropia de buracos negros 
[3,4] e para uma revisão do conceito de entropia termodinâmica de 
buracos negros [5,6]. 


2 * ENTROPIA DE BURACOS NEGROS ATRAVÉS 
DE UMA ABORDAGEM TERMODINÂMICA DE CA- 
MADAS FINAS DE MATERIA 


2.1 « CAMADAS FINAS DE MATÉRIA 


Para se efectuar uma abordagem termodinâmica de um sistema su- 
jeito apenas a acção gravitacional vamos usar camadas finas de ma- 
téria no contexto da relatividade geral que requerem um formalismo 
especial [7,8]. 


e 


O conceito de camada 


na de matéria surge do problema de relativi- 


dade geral sobre como ligar suavemente um espaço-tempo divido em 


duas partes através de uma hipersuperfície. Este problema é resolvido 


deduzindo duas condições que devem ser satisfeitas de modo a que 
o espaço-tempo constitua uma solução válida das equações de Ein- 


stein, algo familiar a outras teorias como o electromagnetismo de 


Maxwell. Estas condições são designadas de condições de junção, e 
estão relacionadas com as propriedades geométricas da hipersuper- 
fície. A primeira delas é 


hal=0 (5) 
onde h,, é a métrica na hipersuperfície - chamada de métrica induzida 
- E os parânteses rectos significam o salto de uma quantidade, isto 
é, a diferença dessa quantidade vista de cada um dos lados da hiper- 
superfície. A primeira condição de junção afirma assim que a métrica 
induzida tem que ser contínua quando se atravessa a ligação das duas 


partes do espaço-tempo. A segunda condição está relacionada com a 


curvatura extrínseca da hipersuperfície, designada de Kab, quantidade 


esta que envolve derivadas segundas da métrica na junção. Se essa 


hipersuperfície infinitesimalmente fina contiver matéria, esta irá pos- 


sui um tensor de energia-momento S.,, digamos. Pode-se mostrar que 


esse tensor e a curvatura extrínseca estão ligados pela equação 


[Kab] — [K]hap = 81Sap (6) 


em que K=hab K.p Esta equação é também conhecida como a equação 
de Lanczos. No caso de não existir nenhuma matéria na hipersuperfi- 
cie, o salto na curvatura extrínseca é nulo, pelo que toda a passagem 
através das duas partes do espaço-tempo se realiza de uma forma 
s 


uave e contínua, como é o caso de uma superfície de uma estrela na 


junção da matéria com o vácuo (para mais detalhes ver [8]). 


Por outro lado, o tensor de energia-momento da camada também 


tem uma forma específica associada à forma como a matéria está dis- 


tribuída. Assumindo que a camada é um fluído perfeito, por exemplo, 


temos automaticamente que S..=(o+p) u, u,+ph,, onde o é a den- 


sidade superficial da matéria na camada, p é a pressão superficial e U 


é a quadrivelocidade da matéria da camada. Através da equação (6), 


dado um espaço-tempo com uma certa métrica, conseguem-se calcu- 


lar as expressões para o e p na camada. 
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2.2 « TERMODINÂMICA DE CAMADAS FINAS 


Suponhamos agora que temos uma camada com uma dada densidade 
superficial o, pressão p e entropia S. Essas quantidades obedecem à 
7º lei da termodinâmica 


TdS = dM + pdA — X;dx; (7) 


onde M é o integral da densidade superficial o - a massa em repouso 
portanto -, À é a área da camada e X, são potenciais generalizados 
correspondentes às eventuais variáveis termodinâmicas x, extras 
que possam existir, como a carga eléctrica por exemplo. A análise ter- 
modinâmica de camada finas no caso não carregado e sem momento 
angular foi feita por Martinez em [9], 


Desta forma, uma simples integração da equação (7) pode levar ao 
cálculo da entropia da camada. Em casos especiais, consegue-se levar 
a camada até ao horizonte de eventos, o que nem sempre é possível, 
Nos casos em que é possível, podemos obter a entropia no limite em 
que levamos a camada para o horizonte de eventos, obtendo dessa 
forma a entropia de um buraco negro. Esta abordagem de levar a ca- 
mada de matéria até ao limite de um buraco negro foi feita pela pri- 
meira em [10,11]. 


2.3 e APLICAÇÃO A CAMADAS 
ELECTRICAMENTE CARREGADAS 


Como exemplo de aplicação deste formalismo, considere-se o caso 
de uma camada esférica electricamente carregada. A métrica origi- 
nada por um corpo de massa m de carga q esfericamente simétrico 
é descrita pelo elemento de linha de Reissner-Nordstróm (ver [1], por 
exemplo) 


=] 
dst=-(1- + Ode +(1-2+0) dr? +r?(do? +sin? 0 dg?) (8) 


que consideramos ser a métrica exterior à camada. No interior da ca- 
mada, a métrica assume-se ser a de um espaço-tempo plano 


ds? = —dt? + dr? + r?(d6? + sin? 0 dg?). (9) 


Assim como no caso da métrica de Schwarzschild se define o raio 
de horizonte de eventos, na métrica (8) surgem dois raios naturais a 


+£ 


definir-se. O primeiro, definido como r,=m+vm2-g2, é o raio gravita- 


cional ou raio de horizonte de eventos da camada carregada, e o se- 
gundo raio é r =m-vVm2-g2, chamado raio de horizonte de Cauchy, que 
se encontra sempre contido dentro do horizonte de eventos. Este úl- 


imo raio é importante para a dinâmica de buracos negros carregados 


mas a sua interpretação física requer uma explicação longa que não 


£ 


será feita aqui, Buracos negros podem ter uma carga mínima nula - 


ecuperando a solução de Schwarzschild com trem - e Uma carga 


máxima q=m, ou equivalentemente r, =r . Neste último caso denomi- 


nam-se de buracos negros extremos, 
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Com as expressões (8) e (9), pode-se calcular as componentes das 
curvaturas extrínsecas e portanto o tensor de energia-momento da 
camada. Considerando que a camada é formada por um fluído perfeito, 
obtêm-se as expressões para a pressão e massa em repouso, que po- 
dem ser substituídas na 1º lei, agora dada por 


TdS = dM + pdA— pda (10) 


em que & é o potencial eléctrico e q é a carga eléctrica. A área é 
A=4mR2, onde R é o raio da camada. 
Da expressão (10) não é Óbvio à partida qual deve ser a forma do po- 
tencial nem da temperatura. No entanto, para o diferencial da entropia 
dS ser exacto, isto é, para se garantir que existe uma expressão para 
a entropia, é necessário que o diferencial (10) satisfaça 3 condições 
de integrabilidade, correspondentes à comutatividade das segundas 


derivadas. Por outras palavras, é necessário obrigar que a derivação 


da entropia em ordem a M e de seguida em ordem a A seja igual a 


derivar em A e depois em ordem a M, por exemplo. Estas relações de 
comutatividade permitem descobrir a forma geral da temperatura e do 


potencial eléctrico, que juntamente com a expressão para a pressão e 
massa em repouso, levam ao resultado 


dS = b(r rn) SEL gr, + blr, 1) Rr gy (11) 


onde b(r,+r Je c(r,tr ) são funções arbitrárias de r,+ e 1. As fun- 
ções b e c são interpretadas como sendo o inverso da temperatura 
da camada e potencial eléctrico, respectivamente, quando esta tem 
um raio infinito. Estas duas funções estão relacionadas por uma outra 
condição de integrabilidade, já que o diferencial (11) também deve ser 
exacto. Note-se, contudo, que este diferencial é independente do raio 
da camada. Por isso, se tomarmos o limite R=r, a forma funcional da 
expressão (11) não muda. 


Para se obter um buraco negro temos que especificar a função b de 


Ee 


uma forma muito particular, dada por 


2 
b(rr) = a (12) 


que é o inverso da temperatura de Hawking. Utilizando esta ex- 
pressão e a condição de integrabilidade de (11), é possível mostrar que 
a solução para a função c é 


c(r,r. ) = a (13) 


Substituindo os resultados (12) e (13) em (11) e integrando, obte- 
mos a entropia 


S= 14) 


em que A, é a área correspondente ao raio gravitacional da camada. 
Dada a escolha (12) e (13), claramente esta entropia é independente 
do raio R da camada. Levando a camada para o horizonte de eventos 
recuperamos assim de uma forma natural a entropia (3) de um buraco 
negro por uma via termodinâmica. Notemos que, após este cálculo, 


verificamos que a entropia de um buraco negro extremo continua a 
ser (3). 


No entanto, pela abordagem termodinâmica conseguimos concluir 
algo mais sobre a entropia de buracos negros extremos. No caso da 
camada, podemos considerar o limite r +=r - antes de se levar para 
o limite de um buraco negro, e só depois calcular a entropia. Supondo 
assim que a matéria da camada satisfaz m=q, o diferencial (11) reduz- 
sea 


dS = f(r.)dr, (15) 


onde f(r, )=b(r. )[-c(r, )r,]. A equação (12) é integrável automati- 
camente sem necessidade de qualquer condição de integrabilidade 
adicional. Por esta razão, f(r,) pode ser escolhida de forma arbitrária, 
já que a função b continua a ter que satisfazer a equação (12) mas 
a função c não está constrangida por nenhuma condição. Conse- 
quentemente, chegamos à conclusão de que a entropia de um buraco 


negro extremo é uma função de r, mas não necessariamente igual a 


(3), a não ser no caso particular em que se escolha a forma (13) para 


c. Este resultado mostra que buracos negros extremos parecem ser 


objectos diferentes de buracos negros usuais, 


3 « CONCLUSÕES 


A entropia de um objecto está relacionada com os graus de liberdade 
microscópicos que o constituem. Num caso de um buraco negro, ex- 
istem dois problemas fundamentais a resolver relativamente à sua 
entropia. O primeiro consiste em saber onde se encontram esses 
graus de liberdade, se no interior do buraco negro ou se no horizonte 
de eventos. O segundo resume-se a averiguar qual a natureza desses 


us de liberdade, o que estes realmente são. Relativamente à pri- 


gra 
meira questão, existe o consenso de que os graus de liberdade se en- 


E + 


contram efectivamente no horizonte de eventos. A nossa abordagem 


ge 


ermodinâmica mostra claramente que estão na camada, uma vez que 


recuperamos a entropia usual de um buraco negro quando a camada 


em um raio igual ao seu raio de horizonte. Quanto ao segundo prob- 


ema, ainda não se sabe a resposta. Em princípio, os graus de liber- 


dade devem ser graus de liberdade gravitacionais, já que um buraco 
negro puro é feito unicamente de gravitação. A nossa abordagem 


ermodinâmica não é suficientemente sofisticada para nos dar essa 


esposta. No entanto, é certamente um problema com raíz na junção 


da mecânica quântica com relatividade geral, pelo que o estudo da en- 


ropia de buracos negros revela ser uma motivação forte na tentativa 


de descobrir uma teoria definitiva de gravitação quântica. 


Para mais detalhes ver [12]. 


| Uma abordagem 
termodinâmica á entropia 
de buracos negros 
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Abstract 


making this issue of extreme interest. 

Other important aspect regarding this theme is that the actual 
national and international bibliography is very incipient and 
limited, regarding the quality of urban runoff in parking lots. 


In a time when water pollution is recognised as an issue of 
major importance and there are just two vears left until the 
deadline imposed by the Water Framework Directive to have 
a good physical, chemical and ecological status in all surface 
water bodies is reached. It rises the importance of studies of 2 * OBJECTIVES 

the potential of urban rainwater to cause significant damages 

inthe receiving water bodies. The objectives developed within this essay were, searching 
This work focuses in the analysis of the rainwater quality in and compiling bibliography and to accomplish a case study 
urban environment, since these waters can drag a significant in Instituto Superior Técnico's exterior parking lot. Diverse 
pollutant load, when coming through the atmosphere and after samples of rain runoff were collected with the purpose of 


flowing over the roofs, roads, and gardens. determinate and evaluate physicochemical parameters. 
Information regarding rain data associated to the events where 
Index Terms sampling occurred was also collected, 


—+ 


— Heavy metals, monitoring, pollution, rainwater, urban runoff All the obtained information was subject to analysis and 
statistical treatment, which allowed the construction of some 
very interesting conclusions. 


1 + INTRODUCTION 
3 « CASE STUDY DEVELOPMENT 


Nowadays no one questions the necessity of the movement 

of the human being. Therefore, independently of the usage, With the purpose of evaluating the quality of urban runoff, 
occasional or professional, motorized vehicles are essential for samplings were made in a specific place of the parking lot. 

the dislocations to be as short and safe as possible. Also, the utilization of certain equipment (udometer), data 
The fact mentioned above makes that, whether you like it or regarding the rain events, such as intensity and duration was 
disapprove it, there are places in cities where naturally parked measured, 
vehicles accumulate, These places can go from simple parking |. The area covered by the study has utilization that involves 
along the streets to parking in proper places such as big silos. directly 53 outside parking places and the access to an 
The waters that result from the washing of the pavement of underground parking area with 101 places. All these are usually 
roads and in specific of parking areas can be, if not properly fully occupied during week time and also by different vehicles 
collected and treated, one of the contamination sources of during the same day. 
superficial water bodies and groundwater, when the discharge. Sampling was made through a simple process, which when 
places of these waters cross with water lines or areas of. the events happened, the water was collected in a plastic 
aquifer recharge, or even through infiltrations that might exist Dox previously cleaned. Then, when the wanted volume was 
a 
S 


ong the path described by the runoff, achieved, the water was transferred to appropriate vessels 
o, these can become a source of pollution, mobile and diffuse, for further analysis in the laboratory. There, the samples were 


mainly of an anthropogenic source. Significant importance is. analysed in the following parameters: 
placed, because of its capability of special dispersion and due 
to the nature of the pollutants in question. “Total suspended solids (TSS); 
Furthermore, in some countries, the diffuse pollution has «Chemical oxygen demand (COD); 
been considered as the main factor able to meet adequate *Biochemical oxygen demand after 5 days 
environmental levels for water quality. incubating at 20ºC (BODS5); 
The fact that parking areas are places where a lot of *Cadmium (Cd); 
manoeuvring happens and where the circulation flows at Copper (Cu); 
very low speeds but sometimes with high revs, powers the “Lead (Pb); 
generation of pollutants, elron (Fe); 

eNickel (Ni); 

eZinc (Zn); 


sTotal hydrocarbons; 
“Oils and fats. 
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At the sampling place, pH, conductivity and temperature were 
measured using appropriate probes. The choice of the param- 
eters for analysis was due to the bibliography analysed. 
Shown below in Table 1 are the characteristics of the shed 
used in this essay. 


Table 1: Values for some parameters of the shed used in the essay. 


Values 


Shed parameters 
Area 230,2 m? 

8,2% in the first 25 meters; 
5,6% in the following 12 meters 
“35 seconds 


Average concentration time 


[3,5-17,4] I/min, according to values of 
precipitation intensity of [0,003-0,014] 
m3/ha.s 


4 « RESULTS AND ANALYSYS 


The following results were obtained. 
ble 2 shows the events where samples were taken and the 
ta associated to them regarding rain measurements. 
able 3a and 3b shows the results returned by the laboratory 
regarding the following parameters: BOD, COD, TSS, cadmium, 
copper, lead, iron, nickel, zinc, oils and fats and hydrocarbons; 
they also show the parameters measures on site: pH, conduc- 
tivity and temperature. 
The analysis of the obtained values (Table 3a and Table 3b) 
allows to verify that the physical parameters remain relatively 
constant, mainly the temperature. 
Very briefly, “sampling” relates to the number of overall sam- 
ples, the “events” relate to the rain events and the relation 
with sampling, and “sample” regards the collected samples per 
event 


Table 2: Description of parameters regarding the events were sam- 
pling was made, 


(mm) 


Precipitation 
duration (min) 
Previous 

dry period (min) 
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Data obtained through the rain bulletin associated to São Julião do Tojal (20C/010) station 
Table 3a: Results from laboratory analysis to the parameters in study. 


Table 3a: Description of parameters regarding the events were 
sampling was made. 
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Table 3b: Results from laboratory analysis to the parameters in 


study. 
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The parameter pH suffered some variations along the pro- 
cess but only between 7,26 (sample 1 from event 1) and 7,98 
sample 1 from event 5). Nevertheless, this small variation is 
extremely relevant, since that values above 7,5, copper, lead, 
inc and cadmium are in less toxic forms for the surrounding 
nvironment, Otherwise, with values bellow 6,5, which were 
ever read during all the study, the toxicity is higher. 
he values of conductivity follow an interesting decrease as 
they allow to verify a phenomena of constant “cleaning” of 
the road, inside events and between events, reaching a value 
of residual pollution around 100 pS/cm after sampling 5.The 
reduction of the conductivity values indicate a possible and 
predictable decrease in the values associated to heavy metals 
along the events. 
It is also important to point out that this “cleaning process” 
was to a certain point visible in the turbidity of the samples 
as along the sampling the water was clearer and clearer after 


collecting each sample. 


— 


DN 


n 
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Regarding TSS, these have an 


after a big previous dry period, returns valu 


and then, it decreases in the 2 
maximum of 3400 mg/l on 


the forth sam 


awkward distributi 
es arou 
samples reaching a 
ple of event 1, de- 


following 


creasing in the following events until the end of t 


It's quite possible that the maxi 


ig discrepancy between val 
tered, 


Despite 


error by the operator during 


the high results of four 


sixth forward, we have again only residual pol 


reduce event after event, with 
annex XVIll from the 


evels below t 


mum of 3400 mg/ 
the sampling proc 
ues and the time | 


th and fifth samp 
ution, which will 


he li 


Portuguese Decree-Law No. 


The values associated with parameters COD and 
not allow the drawing of strong conclusions; the only verifiable 
aspect is that both parameters have a very similar behavior 
within each event. However, the COD values are much higher 
than the BODs for the same samples. 


This fact can be justified by two points: 


on because 
nd 750 mg/l 


he sampling. 
was due to 
ess because 
t was regis- 


es, from the 
mit value of 


236/98º. 
BOD; does 


e Part of the organic matter, difficult to oxidize bio! 
ogically, can be chemically degraded by potassium di 


chromate; 


* Some organic substances can be toxic for the micro 


organisms used 


in the BODs, tests. 


Apart from the above, the fact that there is a garden which 


is regularly taken care of 


response to the high values ins 


ome samples. 


near the sampling place can be a 


Both of the parameters show some samples above the limit 


emission value of annex XVll 
236/982. The ratio between 


of Portuguese Decree-Law No. 
COD and BODs allows under- 


standing the biodegradability of the water. Table 4 presents 


the relation between BODs and 


Table 4 : Relation BOD;/COD. 


COD. 


Event BOD, / COD 
Event | - Sample 2 0,16 
Event | - Sample 3 0,24 
Event 1 - Sample 4 0,15 
Event 5 - Sample 1 0,13 
Event 6 - Sample 1 0,22 
event 8 - Sample 1 0,21 
Event 8 - Sample 3 0,50 
Average 0,26 


Thus, the analysis of the results 
COD shows that only once reg 
representing an event with eas 
were five va 
en0,2 and 
ues below 0,2, indicati 
analysis regarding on 
ts inan average of 0,1 
cative that the waters in 
It 


y the 
Sfort 


is possible to suggest, from 
tially the rai 
ta 

road, there is a decrease in the 


The element zinc, as shown 


ple 5), which presents a slightly 
Iron was the heavy metal that 
tions in all the field work. 
emission level in annex XVII 
first event, and a maximum 


fo 
ine 


Copper and lead showed some 
mit, which is quite positive, as i 
the overall pattern of decrease 
as an additive in fuels. 
A note should also be given to 


ecree-Law no, 236/982, 


Finally, regarding total hydrocar 


concentration range of “accept 


tection limit. 
tis important to point out that 
no. 236/9682 gives respect to 
the ELV shown were the 
ation regarding ELV on 
egislation and European 


Table 6 shows the values fo 
minimum and maximum for the par 


n runoff has higher 
s, and that over the events, with the natural washing of the 


quent increased in biodegradability of the runoff. 


sults were always below the detection limit. 


for the value of the ratio BODs/ 
istered values were above 0,6, 
ily biodegradable waters. There 


ues between 0,2 and 0,6 (4 of them lying betwe- 
0,24), indicating middle biodegradable waters, and 4 
ve of hard biodegradable waters. 


first samples of each event re- 
he relationship BODs/COD, indi- 


itially are hardly biodegradable. 


the presented values that ini- 
concentrations of heavy me- 


inorganic mass, and the conse- 


in Table 3b, has low concentrations 
throughout all the study, with the exception 


of event 4 (sam- 
higher value. 


appeared in higher concentra- 


It showed values above the limit 


rthe first three samples of the 
vent 4 (sampling 5). 


results below the detection li- 
n the case of lead, it's keeping 
since the element was banned 


nickel and cadmium as the r 


For oils and fats, the samples showed very low concentrations, 
| below the limit emission value present at annex XVIII of 


bons, there isn't a limit value of 


discharge that allows the insertion of the results obtained in a 


able”, However, in some events 


it was found that the total hydrocarbons were below the de- 


the annex XVlll of Decree-Law 
residual waters discharge, and 


ones used because there is no legis- 
rainwater discharges on Portuguese 
egislation. 


the average, median, standard deviation, 


ameters. 


Table 6 : Statistical analysis. 


Standard 
Parameters Average Median deviation Min | Max 
15,7 | 22,4 
Tã6 
RE: 
[ET ES REI EM EM 
[07 [05 | 
E 
n 
[Total hydrocarbors mon [ 104 [10] 04 


The result obtained allows the observation of mean values 
higher than the median values. This again shows the variabili 
of the obtained results. This factor is especially evident in 
results obtained for the conductivity, BODS, COD, TSS and iron. 
Regarding the iron and TSS parameters the variability is much 
higher due to the peaks observed in two samples. Therefore, 
it was decided also to show the results for TSS without the 
value that is suspected of being mistaken (3400 mg/). The 
ew results show a much more considerable variability. 

In addition, the values for the standard deviation also give an 
indication of the variability. 


+ E € 
5 
D 


=) 


On Table 7, we have a synthesis of the obtained results or- 
ganized with range and average values and two studies to 
ompare the results, within the comparable parameters, one 
regarding highways and other from urban areas. 


Q 


Table 7 : Final synthesis of the results. 


Concentrations 
Concentrations (mg/1) 
Parameters (mg/1) Urban runoff 
pH 7,61 
(7,26-7,98) 6,5-7,4º 


Concentrations (mg/]) 
Highway A3 Pte Lima? 
ADATS=7849 


Cond (uS/cm) 160,2 
(50,2-362) 33-209 E 
BODs (mg(02)/1) 54,2 
(3-186) 0,3-25 15-141 
COD (mg(02)/1) 299,7 
(6-950) 4,0-256,0 48-964 
TSS (mg/1)º 268,5 
(17-900) 0-159 49-498 
Copper (mg/1) 0,19 
(0,12-0,34) 0-0,045 0,027-0,191 
Iron (mg/1) 2,74 
(0,35-12) 9-3,27 0,071-0,523 
Zinc (mg/1) 0,57 


(0,1-1,3) 0,246-3,839 


0-1 
Oils and fats (mg/1) 2,56 
(0,4-11) 0-52,0 


The obtained results in this study are in the same range as the 
results obtained in the previous presented studies. 

Regarding highways, which naturally have higher traffic and 
speeds quite superior to the ones in parking lots, it is extremely 
interesting to verify that when compared to the values obtai- 
ned in the parking lot the ranges are almost the same which 
indicates that the number of vehicles is mainly independent 
from the results. 


Regarding the values of the urban runoff study, there is not 
much to say because the results are quite similar for the pa- 
rameters where comparison is possible, with the exception of 
iron, which for the parking lot presents values greater than the 
ones on the comparison study. 


Furthermore in Table 8, a correlation table built using the Pe- 
arson coefficient, represents a measurement of the linear 
relation between two quantitative variables. The coefficient 
varies between -1 and 1 where: 


e O means the absence of correlation between variables; 


* -1 means strong linear inverse relation, which means that 
when one variable increases the other decreases and vice- 
-versa; 


e 1 means strong linear relation, which means that when one 
variable increases, the other also increases. 


The parameters subjected to correlation were limited to the 
previous dry period and the parameters analysed in the labo- 
ratory, mainly due to the limited data available, the correlation 
between other parameters didn't return any values of interest 
to the research. 


As the selected parameter was the previous dry period, only 
the first samples of each event were correlated. 

The heavy metals cadmium and nickel didn't report any value. 
Also lead and copper have been removed from the correlation 
table, because they only had values for the first event and one 


value is not enough to correlate data. 


3Values for TSS (mg/]) without the maximum value of 3400 (mg/1) 
Source: Adapted from Barbosa et al, 2011. 

> Average Daily Anual Traffic 
Source: Adapted from Gromaire-Mertz et al, 1999. 


Legend: x - average concentrations; (x-x) concentration range 
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Table 8 : Correlation table between previous dry period and other 
parameters. 


Parameters 


Previous dry period (min) 


COD (me(02)/]) 
TSS (mg/L) 
ron (mg/]) 
Zine (mg/L) 


Total hydrocarbonates (mg/L) 


The analysis of the previous values allowed seeing a strong 
correlation between the previous dry period and the concen- 
tration of oils and fats, BODs, COD and total hydrocarbonates. 
On the other side, a weak correlation was noticed with iron, 
TSS and zinc, 

To conclude, it should be pointed out that no more statistical 
analysis was developed due to the low sampling number, which 
imited the ability to perform deep statistical analysis, thus op- 
ting to do a superficial but secure analysis of the data. 


5 + CONCLUSIONS 


The obtained results in the case study do not differ significan- 
tly from those in other very different situations, corresponding 
o highways with fast traffic. It is also clear that for BODs, COD, 
TSS and iron, the limit emission values constant on the annex 
XVIIl of the Portuguese Decree-Law No. 236/98 are often ex- 
ceeded, 
The results reinforce the idea that this type of water has the 
potential to cause damage to the receiving waters where it is 
ischarged without any treatment, especially when the first 
rainfall events occur after long periods of dry weather; and also 
shows the need to introduce specific legislation for this type 
of water discharge. 
The previous dry period seems to be an extremely influent 
parameter regarding the quality of urban runoff, as this para- 
meter clearly related itself with oil and grease, BODs, COD and 
total hydrocarbons. 
Regarding the heavy metals, nickel and cadmium, they were 
not detected. As for lead and copper, only one value for lead 
and three values for copper were above detection limits. 
tisalso possible to say that the obtained values fall within the 
range of values of most studies associated with urban storm 
water, which presents an organic load usually greater than the 
one recorded in studies associated with highways. 
However, it is interesting to note that the results are within the 
same range of observed values in studies related to highways 
where the number of circulating vehicles is much higher. 
Thus, assuming that vehicles are the main influential factor 
on the quality of urban runoff, the low speeds that are nor- 
mally used inside parking lots seem to be a major factor behind 
the degradation of components of the vehicle and consequen- 
tly the quality of urban runoff. 


+ 


Q 


[04] 


On this point, it seems credible that electric vehicles will bring 
great change, since the maneuvers at low speeds are aided by 
the electric motor, thus avoiding many of the emissions. 
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ERR Modelo Semi-Empírico para 
Vibrações Induzidas por Desprendimento 
de Vórtices em Torres Altas Circulares 


João Parracho 


Resumo 


O objectivo deste trabalho é estudar o problema das vi- 
brações transversais aeroelásticas em torres cilíndricas 
circulares, devidas ao desprendimento alternado de vór- 
tices causado pelo escoamento de um fluido. São obtidas 
amplitudes de deslocamento para as torres usando mé- 
todos analíticos e empíricos. O método principal utiliza- 
do nesta análise é um modelo semi-empírico desenvolvi- 
do no domínio do tempo, proposto por Pinheiro [6]. Este 
método é aplicável ao caso específico de torres isoladas 
de secção circular. Baseia-se em modelos bidimensionais 
onde os efeitos da tridimensionalidade são contabilizados 
através da extensão a condições reais de escoamento e 
da incorporação da análise modal. Estes modelos são o 
modelo empírico não linear de Scanlan [7] e o oscilador 
acoplado de Blevins [3]. Um modelo proposto pela mais 
recente regulamentação europeia (Eurocódigo 1, parte 
1-4 [5]) é também descrito e analisado. Aplica-se o méto- 
do de Pinheiro a casos específicos extraídos da literatura. 
São realizadas comparações com resultados obtidos em 
ensaios de túnel de vento e com outros modelos analíti- 
cos. Estas demonstram uma boa e prometedora concor- 
dância, evidenciando o potencial do método de Pinheiro. 
Abordam-se alguns desenvolvimentos futuros para o mo- 
delo, com o objectivo de estender a sua utilização a um 
âmbito mais alargado de aplicações. Finalmente, apresen- 
tam-se as principais conclusões, juntamente com algumas 
sugestões e precauções para aplicação do modelo. Estas 
demonstram que o modelo de Pinheiro é uma boa ferra- 
menta, sendo melhor usada com cuidado e conjuntamente 
com outras abordagens no projecto de uma estrutura. 
Palavras-Chave: Vibrações induzidas por vórtices; Vibra- 
ções transversais ao vento; Torres circulares; Lock-in; In- 
teracção fluido-estrutura. 


1+ INTRODUÇÃO 


À medida que a tecnologia foi melhorando e materiais 
mais sofisticados foram sendo desenvolvidos, a altura 
das estruturas e a sua esbelteza aumentaram. Estas no- 
vas estruturas têm uma maior flexibilidade e menor amor- 
tecimento, resultado do seu menor peso volúmico. Assim, 
estas estruturas são mais susceptíveis a movimentos in- 
duzidos pelo vento e dimensioná-las para resistir a este 
fenómeno é uma área de crescente importância. As estru- 
turas que exemplificam estes casos são, tipicamente, os 
edifícios altos e as chaminés e torres. 

O vento é a principal acção que induz flexão nestas es- 
truturas, provocando vários efeitos e comportamentos 
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resultantes das vibrações a ele associadas. Devido à 
flexibilidade destas estruturas, podem ocorrer grandes 
amplitudes de deslocamento, velocidades e acelerações, 
podendo provocar danos na estrutura e em partes não es- 
truturais como paredes e fachadas. Em torres metálicas, 
por exemplo, pode ocorrer a redução da sua vida útil de- 
vido à fadiga e os sinais das torres de telecomunicações 
podem ser perturbados devido à vibração das antenas. 
Existem vários fenómenos aeroelásticos, dos quais as ins- 
tabilidades por galope, instabilidade aerodinâmica e as vi- 
brações devidas ao desprendimento alternado de vórtices 
são exemplos. O desprendimento cadenciado de vórtices 
tem sido estudado amplamente ao longo dos anos, sendo 
uma acção susceptível de provocar vibrações de grandes 
amplitudes, principalmente quando ocorre o fenómeno 
chamado de “lock-in”, também conhecido como tranca- 
mento ou captura de vórtices. Este ocorre quando, devido 
à velocidade do vento, a frequência do desprendimento 
é próxima da frequência natural transversal da estrutura, 
originando a separação dos vórtices segundo esta última 
frequência, para uma determinada faixa de velocidade do 
vento. 

Este trabalho foca-se no estudo específico das estruturas 
isoladas de torres e chaminés de secção circular, sob a ac- 
ção de vibrações transversais devidas ao desprendimento 
alternado de vórtices. Para tal, utiliza-se um modelo semi- 
-empírico desenvolvido no domínio do tempo proposto 
por Pinheiro [6]. É essencialmente usado na obtenção de 
estimativas da amplitude transversal máxima esperada 
associada à vibração de torres, quando ocorre o lock-in. 
Na engenharia de vento muitas vezes é suficiente saber 
apenas esse parâmetro para proceder à análise e dimen- 
sionamento da estrutura. 


2 « FORMULAÇÃO DO MODELO 


Ao acontecer o escoamento de um fluido com uma 
determinada velocidade U sobre um corpo não 
aerodinâmico, vai ocorrer uma variação alternada das 
pressões na superfície desse corpo ao longo do tempo. 
Assim, desenvolvem-se forças periódicas sobre o corpo 
quando existe a separação sucessiva de vórtices em cada 
lado do cilindro. Se não existir nenhum obstáculo para o 
movimento do cilindro na direcção transversal à incidência 
do fluido, estas forças vão interagir com a movimentação 
do cilindro, originando as vibrações induzidas por 
desprendimento de vórtices. 

Na dedução do modelo, é assumido que as condições 
de escoamento no cilindro são, aproximadamente, 
bidimensionais. Assume-se que o cilindro tem uma longa 


diâmetro D 


superfície rígida e deslocamentos de corpo inteiro na 
direcção transversal ao escoamento. O cilindro possui 
amortecimento mecânico na mesma direcção, está 
apoiado elasticamente e impedido de se deslocar na 
direcção paralela ao vento. O escoamento é laminar com 
velocidade média uniforme. 


A Figura 1 representa estas condições. 


deslocamento v 


! 


comprimento £ 


Velocidade 
Uniforme U 


Figura 1: Esquema de um cilindro rígido com apoios elásticos 


[B] [6]. 


Se o cilindro também estiver impedido de se deslocar na direc- 
ção transversal ao vento, uma primeira aproximação razoável 
para a força por unidade de comprimento F actuando no cilin- 
dro no tempo t é dada por: 


1 : 
P= 2 PU?DCisen(wrt) (1) 


onde p é a massa específica do ar, U é a velocidade do esco- 
amento, D é a dimensão do cilindro (o diâmetro externo) na 
direcção perpendicular ao escoamento, C, é o coeficiente de 
força lateral e w, = 2mSU/D, sendo S o número de Strouhal. 
No entanto, quando o cilindro pode oscilar na direcção trans- 
versal, esta expressão simples para a força de excitação é 
inadequada. Se y for o deslocamento do cilindro na direcção 
lateral ao vento, então a equação do movimento é escrita na 


forma seguinte: 


my+cy+ky =F(y,9,),t) (2) 


onde m é a massa do cilindro, c a constante de amortecimento 
mecânico e k a rigidez elástica. 

Os modelos apresentados seguidamente são construídos a 
partir da expressão (2), combinando o modelo oscilador aco- 
plado de Blevins com o modelo empírico não linear de Scanlan, 
incorporando depois a análise modal, 


2.1 « MODELO OSCILADOR ACOPLADO DE 
BLEVINS 


Neste modelo, o fluido na esteira é representado como um os- 


cilador não linear auto-excitado, acoplado à estrutura [3]. 

As forças no cilindro são avaliadas de acordo com a teoria da 
dinâmica de fluidos, através da equação fundamental da con- 
tinuidade. Descrições mais detalhadas podem ser encontradas 
na referência [3]. 
As amplitudes de deslocamento máximas na condição de res- 
sonância (lock-in) podem ser expressas pelo produto do rácio 
da massa m com o factor de amortecimento &. Esta variável é 
o amortecimento reduzido 6r, também chamado de número de 
Secruton: 


2m(2m6)  Amng 
s= = oD2 (3) 


As amplitudes de deslocamento máximas para uma estrutura 
cilíndrica no modo natural de vibração | são obtidas, na resso- 


nância, por: 


0,72 
0,3 + RE (| do 
(1,9+6,)S D ( 4) 


Str arms | 


onde y é um parâmetro geométrico que depende da forma mo- 
dal q. A Figura 2 mostra os resultados teóricos deste modelo 
em comparação com resultados experimentais. 
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Figura 2: Amplitudes de deslocamento para cilindros rígidos e 
elásticos em função do amortecimento reduzido [3]. 
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2.2 « MODELO EMPÍRICO NÃO LINEAR DE R. H. 
SCANLAN 


uitas vezes, basta obter as deflexões máximas de uma estru- 
tura. Assim, existe interesse em ter um modelo que reproduza 
alguns dos principais resultados observados perto do tranca- 
mento. 
Para construir tal modelo, é assumido que a excitação, amorte- 
cimento e rigidez aerodinâmicos são fornecidos a um oscilador 
mecânico linear. O modelo adquire a sua não linearidade atra- 
vés da inclusão de um termo cúbico aerodinâmico não linear. 
O seguinte modelo é formulado, em que a equação de movi- 
mento do fluido sob a força induzida por desprendimento de 
vórtices é: 


mlj + 2801) + wty] = 
1 v2y y 
= PU?(D) H(K) (1-cD)p+ dog + (5) 


1 D 
+ 2 CL(K)sen(wt + o) com K = Ea =2ns 


onde H, Ho. C,,£ e a são parâmetros de ajuste obtidos experi- 
mentalmente. H, está relacionado com o amortecimento aero- 
dinâmico linear, H> com a rigidez aerodinâmica e £ é o parâmetro 
não linear do amortecimento aerodinâmico. w, é a frequência 
natural do oscilador mecânico. 
No “lock-in”,C, e H> podem ser desprezados pois são pequenos 
relativamente ao parâmetro de amortecimento aerodinâmico. 
H, e e são determinados a partir da observação experimental 
de duas amplitudes de respostas ressonantes Ay, e Ay, para 
dois valores de amortecimento &, e &>. O valor da amplitude de 


deslocamento Ay é assim: 


1 
RR (pUDH, — 4méw)D]2 
y= pUH,e (5) 


Substituindo Ay,, Ay, € E, E Na equação anterior, temos então 
os parâmetros experimentais, que podem ser usados para pre- 


ver a resposta de uma estrutura em condições semelhantes: 


— BmaS($2 Ay? — 6,472) (7) 
1 D2p(Ay2 — Ay?) 


«= Mi -to)D? 
E Ay — ExAyi 


30 


2.3 « MODELO SEMI-EMPÍRICO DE PINHEIRO 


Este modelo assenta nos dois anteriores, expandindo as suas 
condições bidimensionais para um caso tridimensional, Usa-se 
o conceito de células de vórtices para dividir a torre de secção 
variável sob a acção de vento de velocidade variável U(z), em 
trechos nos quais há desprendimento de vórtices numa deter- 


minada frequência (Figura 3). 


ido :) 


Figura 3: Esquema do modelo analítico de Pinheiro para vi- 
brações induzidas de vórtices em torres esbeltas: (a) células de 
vórtices em torre de secção variável sob acção de vento com 
velocidade U(z), (b) estrutura de secção circular constante sob 
acção de velocidade uniforme (modelo 2D); (c) força aeroelásti- 
ca generalizada F associada ao modo j de vibração [6]. 


Assim, em cada célula | existe uma frequência f, associada a 
valores médios de um diâmetro Di e velocidade Ui relacionados 
com o número de Strouhal da secção de acordo com a seguinte 
equação: 


= (9) 


Através do modelo bidimensional de Scanlan, tem-se que a for- 
ça aeroelástica F(z) em cada trecho é: 


1 ey 9 


+Csen(wet + )] 


onde y é o deslocamento da estrutura em função do tempo. 
Considera-se que cada trecho está inserido numa estrutura 
rígida infinita de secção constante D, sob acção de vento lami- 
nar U, segundo a Figura 3.b. 
Usam-se os resultados teóricos do modelo de Blevins para de- 
terminar as amplitudes de vibração Ay em função do número 
de Scruton, calculando assim H, e e. 
As forças F(z) por unidade de comprimento em cada trecho | 
são então combinadas e transformadas em forças generaliza- 
das, obtendo-se a equação de movimento da estrutura da torre 
em termos da amplitude Y, do seu modo j de vibração: 


Mj(Y, + 2E;w;Y; + wir) =Fm(H,e) (11) 


onde M, &, w; são, respectivamente, a massa, amortecimento 
e frequência natural da estrutura associados ao modo natural 
de vibração j. Fm, é a força modal, função da forma modal q (2) 
e da força por unidade de comprimento F(z). A força modal Fm; 


toma, por sua vez, a seguinte forma: 


Fm; = | eser (ed (12) 


2.3.1 « OBTENÇÃO DOS PARÂMETROS H, E E 


A determinação de H, e £ é feita através de duas amplitudes 
Ay para dois valores de E diferentes, em cada célula com um 
determinado valor de velocidade de vento U,comose demonstra 
na Figura 3, que esquematiza a resposta característica de um 
cilindro rígido de secção constante (modelo bidimensional na 
Figura 3.D). 
Sendo o número de Strouhal para secções circulares em regime 
subcrítico aproximadamente S = 0,2, a velocidade reduzida 
é portanto U,=5(1/S=1/0,2 = 5), valor para o qual são 
determinados, então, os valores máximos de Ay1 e Ay2. 
De acordo com o modelo de Blevins, em cada célula de vórtices 
tiam-se Ay, e Ay, através das equações deste modelo, 
representadas graficamente na Figura 2 em função do número 
de Scruton: 


0,3 + nf 
Ay, = 0,078 — (19 +67) 
Yn =U, (1,9 +6,)S2 cel (1 3) 
AM aiTE 
6, =—n comn=1lor2 
, PDZa 


Considerando os resultados experimentais dados na referên- 
cia [3], admitiram-se respostas lineares nos intervalos 4,0 < 
U,<50e50<U,<6,0, pois através dos modelos estudados 
apenas se conseguem obter os valores máximos. Tal simplifi- 
cação encontra-se representada na Figura 4. Para obter outros 
valores de amplitude para além dos máximos, para velocidades 
reduzidas diferentes de 5, basta fazer uma interpolação linear 
entre os valores já encontrados. 


Ay, é Ay, obtidos 


* da curva do oscilador 
acoplado 


U,=5 


U,- U/fD 


Figura 4: Curvas linearizadas de amplitude em função da ve- 
locidade reduzida para diferentes valores de U, utilizadas na 
determinação dos parâmetros H, e e [6]. 


A faixa de trancamento ou lock-in é definida considerando uma 
variação de 20% em torno da frequência crítica de desprendi- 
mento de vórtices. Fora desta faixa admite-se que o cilindro 
rígido apresenta deslocamentos reduzidos que podem ser 
desprezados. Adopta-se também a simplificação de desprezar 
os resultados para velocidades reduzidas menores que 4 ou 
maiores que 6, uma vez que as amplitudes obtidas são muito 
pequenas e, como tal, irrelevantes para este estudo, apesar de 
na Figura 4 esta consideração não estar totalmente explícita. 
Deste modo, e utilizando os dados anteriores, obtêm-se os pa- 
râmetros H, e £ através das seguintes fórmulas, já escritas em 
cima de outra forma: 


Es Bm ars AyiE, — Ayz& 
Hicet e pD2., ( Ay2 — Ay2 ) (1 4) 
— 4Dia(E, — 2) (15) 


E >>> 
“E AyZE, — Ay2ão 


Seguidamente, é apresentado o caso simples de uma torre de 
secção constante sob vento de velocidade constante ao longo 
da sua altura. 


2.3.2 * TORRE DE SECÇÃO CONSTANTE SOB A 
ACÇÃO DE VENTO DE VELOCIDADE CONSTANTE 
AO LONGO DA ALTURA 


Procede-se à extensão do modelo bidimensional explicado 
em cima, aplicado a uma estrutura de comprimento infinito 
(Figura 5.a), a uma torre encastrada e livre como a da Figura 
5.b, incorporando-se a forma modal. Considera-se, como 
simplificação, que o desprendimento de vórtices se dá numa 
única frequência, o que implica a existência de apenas uma 
célula de vórtices. 


0,2) 


(a) (b) (e) 


Figura 5: Incorporação da forma modal (b) ao modelo rígido (a) 
e correspondente discretização (c) [6]. 
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A força modal pode ser escrita na forma de somatórios modais, 
obtendo-se finalmente: 


1 
Fm, = zPUÉDx 


x E (Domus E emu (16) 


+> Cibjnlesen(wet + a) 
k 


onde Ik são os comprimentos de influência de cada nó k. Segui- 
damente, resolve-se o sistema de equações (11) no domínio 
do tempo pelo método de Runge-Kutta para alguma condição 
inicial. Os deslocamentos em qualquer ponto da torre são obti- 
dos, por sobreposição modal. A correlação espacial também é 
considerada. 


2.3.3 * TORRE DE SECÇÃO VARIÁVEL SOB A 
ACÇÃO DE VENTO DE VELOCIDADE CONSTANTE 
AO LONGO DA ALTURA 


Para O caso mais geral de uma torre de secção variável sob 
a acção de vento de velocidade variável ao longo da altura, 
incorporando a correlação espacial, a expressão da força modal 
tem o seguinte formato: 


1 o 
Fm; = 2º (> O ju Dj HD Us 


-y2 > ge dude) Y 4 (1 7) 


k 


+ a O e CriPjule De UR sen(wit+a) 
K 


onde q é a função de forma da correlação espacial e o subscrito 
i refere-se a cada célula de vórtices. 


3 « DESPRENDIMENTO DE VÓRTICES PELO 
EUROCÓDIGO 


O Eurocódigo fornece um método simples para estimar as vi- 
brações induzidas por desprendimento de vórtices em estrutu- 
ras. Neste trabalho, opta-se apenas pela aplicação do Método 
1, uma vez que é um método mais geral que pode ser utilizado 
para vários tipos de estruturas e de configurações modais. O 
deslocamento máximo, Ye max é calculado pela seguinte expres- 
são; 


Yemax 1. 1 
b St? Sc 


ES Kw” Ciat (18) 


Se 


onde b é o diâmetro exterior do cilindro de base circular, St O 
número de Strouhal, Sc o número de Scruton, KW o coeficien- 
te de comprimento de correlação efectivo, K o coeficiente de 
configuração modal e clat é o coeficiente de força lateral. Expli- 
cações detalhadas podem ser encontradas no EC 1-4, anexo E, 
Nota-se que a metodologia proposta pelo Eurocódigo é sim- 
plista, não considerando vários parâmetros que podem in- 
fluenciar a resposta e podendo originar resultados distantes 
da realidade. As amplitudes deste método são em geral baixas 
comparativamente a outros métodos. Deve ser então aplicado 
com uma análise cuidada e após proceder à sua correcta cali- 
bragem, particularmente em relação à taxa de amortecimento. 


4 « RESULTADOS 


4.1 « DESCRIÇÃO GERAL DOS EXEMPLOS 


As estruturas examinadas seguidamente são, tal como refe- 
rido, apenas estruturas cilíndricas com secção circular, pois é 
esse o domínio de aplicação do modelo estudado. São anali- 
sadas uma haste cilíndrica com menos de 1 m de altura e duas 
chaminés com alturas a rondar os 100 m. As estruturas são 
compostas por materiais diferentes, variando também o nível 
de amortecimento, os perfis de velocidade do vento aplicados 
e as dimensões das secções transversais. 


4.2 * EXEMPLO 1 - HASTE DE MADEIRA EM 
TUNEL DE VENTO 


Neste exemplo procede-se a uma comparação entre os 
resultados obtidos em túnel de vento para uma haste 
cilíndrica de secção constante, segundo a referência [2], e os 
resultados teóricos calculados pelos métodos de Pinheiro e do 
EC 1-4. A Figura 6 representa a esquematização do modelo 
de teste utilizado no túnel de vento e as dimensões da secção 
transversal da haste cilíndrica. 
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Figura 6: Esquema da haste e respectivas dimensões [2]. 


A haste é feita de madeira de densidade pmadeira + 240 kg/m'. 
A massa por unidade de comprimento é m = 0,77 kg/m. O seu 
topo está rotulado e a rigidez foi aumentada acoplando trans- 
versalmente duas molas lineares na sua base livre, cada uma 
com rigidez K = 1968 N/m. A frequência natural é f = 20,2 Hz 
e a taxa de amortecimento é E = 0,0031. Foram aplicados um 
perfil de velocidade constante e outro linear. O perfil constante 
foi aplicado com incrementos sucessivos de velocidade e o per- 
fil linear imposto com a velocidade máxima a ocorrer na extre- 
midade livre inferior do cilindro, procedendo a incrementos ou 
decrementos progressivos da velocidade. O perfil da velocidade 
variável é dado pela expressão (19). Os resultados obtidos são 
apresentados na Figura / e Figura 8. 


V(z) = Vo(1 + 0.887) (19) 


onde V é a velocidade em z = 0, 
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Figura 7: Amplitude de resposta do cilindro sob perfil cons- 
tante de velocidade do vento, para o primeiro modo de vibração 
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Figura 8: Amplitude de resposta do cilindro sob perfil linear de 
velocidade do vento, para o primeiro modo de vibração. 


Para vento constante, as amplitudes adimensionais máximas 
obtidas no modelo de Pinheiro são de 0,44 e 0,25 para correla- 
ção total e parcial, respectivamente, ambas para U, = 5. Compa- 
rando com o túnel de vento, verifica-se uma boa conformidade 
ntre ambos, com ligeiras diferenças nas amplitudes e velo- 
idades reduzidas para as quais aquelas ocorrem. A faixa de 
rancamento de Pinheiro é mais larga que a experimental. Estas 
iferenças devem-se a alguns factores como disparidades no 
úmero de Strouhal, com a consequente influência na velocida- 
e crítica e nas frequências naturais, e até devido à discretiza- 
do da haste no programa de elementos finitos, influenciando 
orma modal obtida. A correlação total é um limite superior 
posta, com melhores resultados neste caso específico, 
uma hipótese mais fidedigna. A máxima amplitude ex- 
ntal corresponde à média das amplitudes máximas para 
as correlações. 
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A amplitude máxima segundo o Eurocódigo é 0,46, quase igual 
à do método de Pinheiro para correlação total, embora ocorra 
para uma velocidade reduzida maior. A forma do gráfico é dife- 
rente dos outros porque a formulação do EC 1-4 consiste numa 
simples multiplicação de parâmetros. Estes aumentam com o 
crescimento da velocidade do vento até um certo limite a partir 
do qual permanecem fixos, e daí a amplitude ser igual para ve- 
locidades reduzidas maiores que 6,7. 

Para vento linear, as amplitudes máximas são 0,43 e 0,25 para 
correlação total e parcial, respectivamente, para U, = 5,70. 
Existe concordância razoável nos resultados, sobretudo na 
correlação total, pois é uma hipótese mais adequada devido 
à bidimensionalidade do escoamento, exemplificado pela sec- 
ção transversal constante. Não existem resultados do EC 1-4 
para velocidade linear porque o Método 1 implica que o cilindro 
esteja sujeito a velocidade constante, pois apenas se querem 
calcular os deslocamentos máximos. 


4.3 + EXEMPLO 2 - CHAMINÉ DE BETÃO 
ARMADO EM ITÁLIA 


Comparam-se os resultados obtidos com um método numérico 
referenciado em [4] com os calculados pelos modelos de 
Pinheiro e do EC 1-4, A chaminé em estudo está representada 
na Figura 9. 


Figura 9: Chaminé de betão armado em Trieste, Itália [4] 
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A estrutura de 100 m de altura é cilíndrica, com secção trans- 
versal constante de diâmetros exterior e interior de 63 me 
5,7 m, respectivamente, A densidade é pperso = 2500 kg/m? e 
a massa por unidade de comprimento m = 14137 kg/m. A fre- 
quência natural é f = 0,37 Hz. Foi aplicado um perfil constante 
de velocidade ao longo da altura, sendo aumentado progressi- 
vamente à semelhança do exemplo anterior, variando apenas 
o coeficiente de amortecimento & Apenas foi considerada a 
correlação total, pois o escoamento é fortemente influenciado 
por condições bidimensionais. Os resultados são apresentados 
nas figures seguintes. 
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Figura 10: Amplitude de resposta da chaminé sob velocidade 
constante para o primeiro modo de vibração, com E = 0,01. 
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Figura 11: Amplitude de resposta da chaminé sob velocidade 
constante para o primeiro modo de vibração, com & = 0,005. 
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Figura 12: Amplitude de resposta da chaminé sob velocidade 
constante para o primeiro modo de vibração, com E = 0,01. 
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O intervalo de trancamento da Figura 4 foi alterado para estar 
melhor ajustado aos resultados numéricos. Considera-se que o 
trancamento se inicia em U, = 4 e termina em U, = 7, como pico 
ocorrendo em U, = 5,5, 


As amplitudes adimensionais máximas para o modelo de Pi- 
nheiro são de 0,0518, 0,1092 e 0,5319, para E = 0,01, E = 
0,005 e & = 0,001, respectivamente, ocorrendo para U, = 5,5. 
Existe boa concordância entre os resultados. As discrepâncias 
observadas são explicadas pelas mesmas razões do exemplo 1. 
Verifica-se que, conforme o amortecimento diminui, as diferen- 
ças entre as amplitudes dos dois métodos vão aumentando. Tal 
acontece porque o método numérico tem maiores capacidades 
auto-limitativas da amplitude em condições críticas não amor- 
tecidas, comparativamente ao método de Pinheiro. A forma do 
gráfico do EC 1-4 é equivalente ao do exemplo 1, mas os resul- 
tados são bastante diferentes, sendo cerca de um terço dos 
outros dois métodos. Este método é muito sensível à grandeza 
do amortecimento & podendo dar resultados bastante mais 
baixos. Conclui-se que o EC 1-4 deve geralmente aplicar-se 
com valores mais conservativos de amortecimento, requerendo 
um ajuste e calibragem cuidadosos ao efectuar a sua compara- 
ção com outros métodos. 


4.4 + EXEMPLO 3 - CHAMINÉ DE AÇO NA 
TAILÂNDIA 


Comparam-se os resultados obtidos com o modelo de Vickery- 
Basu referenciado em [1] com os calculados pelos modelos de 
Pinheiro e do EC 1-4, A chaminé analisada está representada 
na Figura 13. 
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Figura 13: (a) Chaminé de aço com 90 m de altura em 
Rayong, Tailândia. (b) Dimensões básicas da chaminé [1]. 


A chaminé tem 90 m de altura e secção transversal cilíndrica 
sem revestimentos. É composta por quatro secções soldadas 
entre si, com diâmetros e espessuras diferentes, conforme a 
Figura 13.b. O aço tem densidade paço = /850 kg/m3, e a massa 
média por unidade de comprimento é m = 1907 kg/m, sendo 
a massa total de cerca de 177000 kg. A frequência natural é 
f = 0.93 Hz, Como as secções constituintes da chaminé têm 
diâmetros diferentes, logo têm velocidades críticas diferentes, 
Para obter então a máxima resposta de amplitude, foram apli- 
cadas velocidades de vento diferentes nas secções, mas cons- 
tantes ao longo de cada uma. Assim, para U, = 5 em todas as 


secções, foram aplicadas velocidades de 2418 m/s, 19,07 m/s, 
14,88 m/s e 10,23 m/s para as secções desde a base até ao 
topo, respectivamente. 


Apenas as amplitudes máximas são comparadas, pois são os 
resultados mais importantes e em [1] apenas se indicam estes 
valores, sem qualquer menção às velocidades reduzidas para os 
quais são obtidos. Os resultados são apresentados na tabela 


seguinte. 
L] 


Amort. (€) 


Amplit. máx. de deslocamento (m) 


Pinheiro Vickery-Basu 


039 


TABELA 1: (Amplitudes absolutas máximas de deslocamento 
da chaminé, para várias taxas de amortecimento. 


As amplitudes têm boa concordância, excepto para o método 
do EC 1-4. Novamente, à medida que o amortecimento dimi- 
nuli, maiores são as amplitudes. Este crescimento da amplitude 
é exponencial, por razões já mencionadas anteriormente. Os 
resultados do EC 1-4 são similares aos exemplos precedentes, 
sendo cerca de um décimo dos calculados pelos outros dois 
métodos. A razão para tal acontecer já foi explicada. 


5 e SUGESTÕES PARA TRABALHOS FUTUROS 


o) 


O modelo de Pinheiro demonstra ter potencial para ser aplicad 
a situações mais gerais e complexas. Estas podem passar 
pela consideração de outros regimes de escoamento para 
além do subcrítico, tais como o crítico e o supercrítico, ou pe 
implementação de diferentes configurações realistas de célula 
de vórtices. A extensão do modelo a outros tipos de estrutura 
e secções transversais, para além das torres e chaminés e d 
secção circular, é igualmente interessante. 


fab) 


ODE VA RA) 


6 « CONCLUSÕES 


Os efeitos do vento perpendicularmente à direcção da sua 
actuação nas estruturas é um assunto complexo que ainda 
está longe de ser totalmente compreendido. 
Foi apresentado um modelo semi-empírico desenvolvido por 
Pinheiro para a estimativa de amplitudes de deslocamento 
transversal em torres de secção circular sob a acção do vento, 
quando ocorre o fenómeno de desprendimento de vórtices. 
Foram comparados resultados experimentais de modelos 
reduzidos ensaiados sob escoamento laminar em túnel de 
vento e resultados obtidos por outros métodos numéricos e 
teóricos. Nos exemplos estudados verificou-se existir uma boa 
concordância entre todos os resultados, apesar das dificuldades 
e incertezas que existem na análise dos problemas. Apenas 
em alguns casos houve a necessidade de proceder a uma 
calibragem e ajuste dos parâmetros de entrada do modelo de 
Pinheiro, como, por exemplo, através da alteração da janela de 
lock-in predeterminada. 


Já a comparação com o método do Eurocódigo produziu 
resultados insatisfatórios. Apenas no exemplo 1 as amplitudes 
máximas obtidas pelos dois métodos foram praticamente 
iguais, o que demonstra que é um modelo que deve ser utilizado 
com precaução, sem descurar outros métodos que possam 
eventualmente ser mais exactos. Tal como referido, a escolha 
do amortecimento deve ser muito criteriosa e ponderada, pois 
é um parâmetro de extrema relevância no EC 1-4. 


Conclui-se que o modelo de Pinheiro constitui uma ferramenta 
bastante útil na análise do fenómeno de desprendimento 

de vórtices em torres circulares. O modelo foi utilizado 

numa forma elementar mas tem simultaneamente bastante 
potencial para ser estendido a casos mais gerais e amplos, 
desde que se conheçam determinados parâmetros 
fundamentais. 


Como tal, o método estudado neste trabalho é um 
instrumento bastante vantajoso, mas que não deve ser nunca 
utilizado isoladamente. Deve ser antes utilizado em conjunto e 
em complementaridade com outras ferramentas em diferentes 
fases do projecto, de modo a assegurar uma análise exaustiva, 
fiável e de qualidade. 
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Leonor Guedes da Silva 
A y À 


Estudante do mestrado em Engenharia Biológica do IST, Leonor Guedes 

da Silva coordenou uma equipa de investigação que venceu a edição 2013 

do Prémio do Melhor Projecto do Mestrado de Tecnologia e Ciências da Vida da Delft University of Technology. 
O projecto baseava-se no processo bioquímico que permite alterar os actuais processos de produção de energia por queima directa 
de pellets de madeira, transformando-os em etanol e ácido succínico, reduzindo em cerca de 44% as emissões de C02 para a 
atmosfera. A Técnica falou com a Leonor sobre a experiência e de como esta se irá repercutir no seu futuro profissional. 


Porquê Engenharia Biológica? Sempre foi clara a ideia de seguir a área das 
ciências? 


Apesar do meu fascínio pelo ser humano (que naturalmente me atraía para Medicina), o Instituto Superior 
Técnico sempre foi a minha primeira escolha de escola por ser uma instituição de excelência com várias 
ligações a universidades estrangeiras que fazem com que este seja, no fundo, “uma ponte para o mundo”, 


A grande questão que então se pôs foi: Como conciliar este fascínio e a vontade de pertencer ao 
Técnico? Felizmente, encontrei a minha resposta em Engenharia Biológica; um curso versátil, que 
combina a matemática, a química, a física e a biologia num só e que me permitirá aceder a um vasto leque 
de saídas profissionais desde a investigação à industria, passando também por gestão e consultoria. 


Porquê o Instituto Superior Técnico na tua escolha e em que medida 
consideras que foi uma mais-valia em termos curriculares? 


Enquanto estive no estrangeiro, apercebi-me que os alunos do Técnico possuem tantas ou mais 
capacidades quanto alunos vindos de universidades do topo dos rankings mundiais. Esta percepção 

é bastante reconfortante e penso que este facto se deve à excelente “bagagem” académica e 
extracurricular que o Técnico providencia aos seus alunos e que, apesar da atual situação portuguesa, é 
o principal motivo que mantém elevada a empregabilidade dos alunos diplomados desta instituição. 
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Leonor Guedes da Silva 


Como escolheste o país para Erasmus? Ou seja, quais os 
aspectos que te levaram a tomar a tua decisão? Pensaste 
antecipadamente nas novas oportunidades que lá poderias 


encontrar? 


Dado que fui em 
entre o meu cur 
programa de Eng 


com óptimas con 


Erasmus um an 
so e a universi 


dições e que, e 


enharia Bioquímica sem 
estava garantida. Para além deste facto 
m conjunto com a industria, faz 
na minha área. Para além disso, está situ 
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r, TU Delft é uma universidade de excelência, 


ada na Holanda, Um país onde a q 


investigação à indústria, passando ta 


mbém por gestão e consultoria. 


bastante investigação 


ualidade de 
ais desde a 


Como foi a estadia em Delft em termos académicos? 


Em termos académicos, TU Delft excedeu as minhas expectativas. Para além do que 
mencionei anteriormente, permitiu-me ver o conceito “universidade” de uma maneira 
diferente, onde o balanço entre as atividades académicas e extracurriculares não é 
só importante, como é encorajado, dado que para o empregador é essencial que o 
candidato possua não só competências técnicas mas como também soft-skills (saber 
comunicar, trabalhar em grupo, etc,). 


Fala-nos um pouco do projecto que integraste enquanto 
estiveste a fazer Erasmus, de que se tratou esta 
investigação? Como surgiu esta ideia inovadora? 


O objectivo do meu professor era juntar grupos de alunos, com diversos backgrounds 


dentro da bioen 


de negócios apo 
integrada de bi 
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ocombustíveis-bioquim 
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genharia, para trabalharem juntos num caso real, relevante para a 
nternacional, transformando o seu conhecimento num plano 
iado por um projeto técnico-económico conceptual para 
icos-bioenergia (neste caso, produção 
e captura de CO, para produção de químicos 
rescentado). No final, os planos de negócios foram apresentados 
a um júri de indústrias, professores e investidores. O nosso foi o melho 


a produção 
de 


re por isso 


Como foi trabalhar em equipa e coordenar um projecto tão 
vanguardista? 


Foi sobretudo desafiante, pois para além de termos de assumir responsabilidades 
perante terceiros (cliente), tínhamos a motivação adicional de conseguir uma 
boa classificação numa das unidades curriculares mais importantes do mestrado. 
A tarefa que me coube foi muito facilitada pela excelente qualidade do grupo 
(que resultou de um sorteio) e pelo elevado empenhamento de todos, sobretudo 
cumprindo o plano que elaborámos para a distribuição de tarefas. 


Quais os aspectos que te fascinaram 
durante esta investigação e, pelo 
contrário, quais os maiores desafios que 
tiveste de superar? 


Para mim, este foi um projeto aliciante por tratar de um 
problema de grande atualidade e motivador por sentir 
que estava a contribuir para a resolução de um problem 
que nos afecta a todos (sobre-exploração de recursos nã 
renováveis e consequente emissão excessiva de (02). Par 
além disso, é sempre um desafio tentar resolver problema 
que ainda ninguém resolveu, mas em contrapartida 
dificuldade é acrescida dada à falta de publicações sobr 
o tema. 
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Entretanto, já recebeste propostas de 
financiamento para dar continuidade 
ao projecto? Gostarias de continuar a 
desenvolver e aprofundar formas de 
substituir os petroquímicos através da 
utilização de outros resíduos? 


O design ainda está apenas na fase conceptual pelo qu 
terá de ser aprofundado por engenheiros e/ou alunos d 
doutoramento profissional em engenharia (um program 
holandês para preparação de mestres para o trabalh 
em industria baseado em design aplicado à engenharia 
Dentro ou não deste projeto, gostaria bastante d 
continuar a contribuir para o desenvolvimento de um 
“bio-economia” sustentável que se baseia na produção 
de químicos, materiais e energia a partir de bio-recursos 
renováveis em vez do tradicional petróleo. 
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Depois de venceres este prémio 
surgiram-te novas oportunidades? 
Propostas de integração noutras 
investigações? 


A noção de oportunidade depende da perspectiva de cada 
um de nós; no meu caso, não tive propriamente ofertas 
de emprego mas tenho a certeza de que o conhecimento 
adquirido durante o projeto e as ligações estabelecidas 
abrir-me-ão novas oportunidades às quais, possivelmente, 
sem ele não teria acesso, 


Projectos futuros, já tens algo em 
mente? Pretendes seguir uma carreira na 
área da investigação científica? 


Ainda tenho quase um ano pela frente até acabar o meu 
mestrado, por isso muita coisa pode mudar. Gostaria de 
fazer uma breve passagem pela industria para sair um 
pouco do meio académico e “ganhar outra perspectiva” e 
depois gostaria de fazer um doutoramento. 


Do teu ponto de vista, qual/quais 0/s 
aspecto/s mais recompensadores de 
trabalhar no âmbito da investigação? 


Trabalhar em investigação não é fácil, não é como fazer um 
exame em que a matéria vem nos livros, é preciso ser-se 
criativo, paciente e persistente. A recompensa do trabalho 
árduo de um investigador aparece com as “pequenas 
conquistas” feitas em direção a uma possível solução do 
problema em estudo, especialmente quando o problema 
tem um grande impacto na sociedade. 


Quais os conselhos que darias aos 
estudantes do IST que gostariam de 
envergar pela área da investigação? 


Antes de decidirem se gostariam ou não de seguir pela 
investigação, ou até mesmo que área de investigação 
gostariam de seguir, experimentem investigação em 
diversos temas, aproveitem as oportunidades criadas pelos 
vossos núcleos de alunos (como por exemplo o programa 
de estágios do NEBIST) e alarguem os vossos horizontes. É 
importante descobrirem o que vos realizará mais na vossa 
vida porque “quem corre por gosto, ndo cansa”, 
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O Centro de Estudos em Gestão do Instituto Superior 
Técnico (CEG-IST) é uma unidade de investigação em Gestão 
ntegrada numa escola de Engenharia, dedicada à análise de temáticas 
que promovem a interacção entre gestão e tecnologia. A investigação 
em sistemas de engenharia e gestão procura analisar problemas 
relacionados com a gestão de sistemas e tecnologias utilizando 
métodos quantitativos para formular problemas complexos e propor 
soluções robustas que sejam aplicáveis a decisões e estratégias 
empresariais e de políticas públicas. 


A missão do CEG-IST consiste em: 
- Criar, desenvolver e partilhar novo conhecimento; 


- Validar o conhecimento criado junto da comunidade 
académica e da sociedade em geral; e 


- Contribuir para a formação de investigadores 
e para a melhoria do ensino universitário. 


Estabelecido em 1999, o CEG-IST integra actualmente 31 
investigadores doutorados e mais de 40 alunos de investigação. Os 
investigadores do Centro focam primordialmente as seguintes áreas: 
- Análise de Decisão; 


- Operações e Logística; 


- Optimização; 


- Análise de Desempenho; e 


- Inovação, Empreendedorismo e Desenvolvimento. 


O CEG-IST desenvolveu um sistema interno de avaliação da actividade 
entífica dos seus investigadores baseado na sua produtividade e na 
ualidade das publicações. Este sistema serve de guia para a afectação 
e recursos financeiros às actividades científicas dos membros. As 
blicações recentes dos membros do CEG-IST incluem artigos numa 
versidade de revistas internacionais, incluindo, por exemplo: Journal 
fthe Operational Research Society; Omega; Energy Policy; Journal of 
roductivity Analysis, Research Policy, European Economic Review; and 
IT Evolutionary Computing Journal. 
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O Instituto: 

Criado no âmbito do Programa CIÊNCIA em 1992, o IDMEC - Instituto de Engenharia Mecânica é uma 
Associação privada sem fins lucrativos, tendo como membros fundadores o Instituto Superior Técnico (IST) 
da Universidade Técnica de Lisboa, e a Faculdade de Engenharia da Universidade do Porto (FEUP). 


O IDMEC é constituído por dois Pólos: em Lisboa, no Instituto Superior Técnico, e no Porto, na Faculdade 
de Engenharia da Universidade do Porto, dotados de autonomia administrativa e financeira, sendo que cada 
Pólo tem órgãos científicos e de gestão próprios. Cada Pólo está organizado em Centros de Investigação: 


Pólo do IST 

Centro de Projecto Mecânico 

Centro de Sistemas Inteligentes 

Centro de Tecnologia Mecânica e Gestão Industrial 
Centro de Energia e Mecânica de Fluidos 


Em Junho de 2008 foi criado o LAETA - Laboratório Associado de Energia, Transportes e Aeronáutica. 
O LAETA é um consórcio liderado pelo IDMEC, com 6 Unidades de Investigação e 230 Doutorados, sendo 
que 98 são do IDMEC- Pólo IST. 


Unidades de Investigação do LAETA: 

e Instituto de Engenharia Mecânica, Pólo do Instituto Superior Técnico 

e Instituto de Engenharia Mecânica, Pólo da Faculdade de Engenharia da Universidade do Porto 

e Instituto de Engenharia Mecânica e Gestão Industrial, Faculdade de Engenharia da Universidade do Porto 
e Laboratório de Aerodinâmica Industrial, Faculdade de Ciências e Tecnologia da Universidade de Coimbra 

e Centro de Ciências e Tecnologias Aeronáuticas e Espaciais, Instituto Superior Técnico 

e Centro de Aeronáutica e Astronáutica, Universidade da Beira Interior 


Objectivos e Missão 

Os objectivos do IDMEC são o exercício de actividades de investigação científica fundamental e aplicada, de 
desenvolvimento experimental, de inovação tecnológica, de formação profissional, de pós-graduação e de 
prestação de serviços no âmbito da engenharia mecânica. 


A missão do LAETA desenvolve-se em torno da investigação científica e tecnológica em consórcios de 
investigação nacionais e internacionais, da transferência de tecnologia para as empresas e da formação de 
engenheiros altamente qualificados e com uma forte exposição a I&D de base tecnológica. 


Governação: O IDMEC/IST é gerido por uma Comissão Coordenadora, composta pelo Presidente e Vice- 
Presidente do Instituto e pelo Coordenador de cada Centro; 


Prof. Carlos A. Mota Soares > (Presidente) 

Prof. José Carlos Pereira > (Vice Presidente) 

Prof. Hélder C. Rodrigues > (Coordenador do Centro de Projecto Mecânico) 

Prof. Paulo Firme Martins > (Coordenador do Centro de Tecnologia Mecânica e Gestão Industrial) 
Prof. José Sá da Costa > (Coordenador do Centro de Sistemas Inteligentes) 

Prof. Mário Costa > (Coordenador do Centro de Energia e Mecânica de Fluidos) 
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O Instituto de Sistemas e Robótica (ISR-Lisboa) é uma unidade 
de investigação e desenvolvimento (l&D) de base universitária onde 
sdo desenvolvidas atividades de investigação multidisciplinar avançada 
nas áreas da Robótica e Processamento de Informação, incluindo Teoria 
de Sistemas e Controlo, Processamento de Sinal e Imagem, Visão por 
Computador, Sistemas Biomédicos, Inteligência Artificial e Sistemas 
Inteligentes. Resultados dessa investigação são aplicados em áreas 
como os sistemas robóticos aéreos, terrestres e oceanográficos 
autónomos, robôs de busca e salvamento, robôs auxiliares que 
interagem naturalmente com humanos, comunicações móveis, sistemas 
multimédia e imagem médica. 


Ex-estudantes e investigadores do ISR-Lisboa são atualmente 
professores em universidades norte-americanas e europeias, ou têm 
posições seniores em empresas nacionais e internacionais, tendo ainda 
criado 6 spinoffs (IdMind, Observit, Blue Edge, Albatroz, Self Tech, 
Reverse Engineering) em áreas afins às da sua I&D, 


O ISR-Lisboa tem um registo antigo e notável de atividades de 
divulgação e promoção da cultura científica, nomeadamente usando a 
robótica para incentivar os jovens a estudar ciências e tecnologia e a 
seguir carreiras científicas, bem como para demonstrar a importância 
da investigação na educação para a vida, tendo estado na vanguarda 
de uma infinidade de projetos e atividades do Ciência Viva, incluindo 
ocupação científica de jovens nas férias e projetos com escolas do 
ensino secundário e profissional. 


Desde 1996, o ISR-Lisboa é sujeito a avaliações periódicas, por 
painéis internacionais de cientistas independentes organizados pelo 
Ministério da Ciência, Tecnologia e Ensino Superior, tendo sempre 
obtido a classificação máxima de Excelente. Em conjunto com o IN 
+,0 CAMGSD, unidades de I&D do IST, o CREMINER da Faculdade de 
Ciências da Universidade de Lisboa, o Centro do IMAR da Universidade 
dos Açores e o MITI da Universidade da Madeira, o ISR-Lisboa faz parte 
desde 2011 do Laboratório Associado LARSyS (Laboratório de Robótica 
e Sistemas em Engenharia e Ciência). 
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Crédito Formação BPI/AEIST 


O seu parceiro 
de estudo, ||| 


Invista na sua formação com o BPI. Com prestações mensais muito reduzid 
o Crédito Formação BPI/AEIST (TAEG de 4,2%)! oferece condições especia 
aos sócios da AEIST, para financiar as suas despesas com Licenciaturas, 
Pós-Graduações, Mestrados ou Doutoramentos, em Portugal e no estrangeiro, 
nomeadamente pagamento de propinas, material didáctico, viagens e estadias | 


mento de € 30.000 
idexada a Euribor 


ITAEG de 4,2%, mensalidade de € 302,67 para um fil 
a 120 meses, sem período de carência e com taxa 
a 3 meses (0, 274%). pe de 3,5%. Montante totz 
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